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Abstract 

This is the first of a series of papers devoted to the initial value problem for the Euler system of 
compressible fluids and augmented versions containing higher-order terms. We encompass solutions 
that have finite total energy and enjoy a certain symmetry (for instance, plane symmetry); these solutions 
may have unbounded amplitude and contain cavitation regions in which the mass density vanishes. In 
the present paper, we are interested in dispersive shock waves and analyze the zero viscosity-capillarity 
limit associated with the Navier-Stokes-Korteweg system. Specifically, we establish the existence of finite 
energy solutions as well as their convergence toward entropy solutions to the Euler system. We encompass 
a broad class of nonlinear Navier-Stokes-Korteweg constitutive laws, which is determined by two main 
conditions relating the viscosity and capillarity coefficients, that is, on one hand the strong coercivity 
condition (as we call it) which provides a favorable sign for the integrated dissipation associated with 
an effective energy, and on the other hand the tame capillarity condition (as we call it), which restricts 
pointwise the strength of the capillarity relatively to the viscosity. Rather mild conditions on the growth 
of the constitutive functions are aso imposed, which are required in order to define finite energy weak 
solutions to the Navier-Stokes-Korteweg system, even in the presence of cavitation. Our method of 
proof relies on fine algebraic properties of the Euler system and combines together energy and effective 
energy estimates, dissipation and effective dissipation estimates, a nonlinear Sobolev inequality, high- 
integrability properties for the mass density and for the velocity, and compactness properties based on 
entropies. 
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1 Introduction and outline 

1.1 Euler system of compressible fluid flows 

In this paper and the companions | |24l [25)1 , we establish existence, compactness, and convergence results 
for (augmented versions of) the Euler system of compressible fluid dynamics, when the fluid flow under 
consideration has plane, cylindrical, or spherical symmetry or the fluid evolves within a nozzle with 
variable cross-section. The proposed Finite Energy Method, as we call it, encompasses real fluids (as the 
standard class of polytropic fluids is too restrictive in applications) as well as broad classes of augmented 
Euler systems, which incorporate physically-relevant small-scale terms. The proposed method allows 
us to validate the singular limit problem associated with these models — in presence of cavitation and 
shock waves — and, specifically, to establish the convergence with finite energy solutions to the augmented 
models toward finite energy solutions to the Euler system. It originates from pioneering works by DiPerna 
||T81[19 | on bounded solutions and by LeFloch and Westdickenberg [47J on finite energy solutions. 

In the present paper, we rigorously validate the vanishing viscosity-capillarity method under mild and 
physically realistic assumptions, as is now presented. Roughly speaking, we determine conditions under 
which 'dispersive shock waves' converge to shock wave solutions to the Euler system when the capillarity 
and viscosity tend to zero. Recall that the Euler system for isentropic compressible fluid flows (in plane 
symmetry, for simplicity) reads 

Pt + (P").v = 0' 
{pu)t + ipu^ + pip)), = 0, 

in which t > represents the time variable while the spatial variable x takes its values in either the torus T 
or the real line ]R. The fluid is characterized by its mass density p = p{t, x) > and velocity u = u{t, x) e R. 
The two equations l Il.Hl form a nonlinear hyperbolic system of conservation laws, provided the pressure 
function p = p{p) satisfies the monotonicity condition 

p'ip) > (p > 0). (1.2) 

Importantly, strict hyperbolicity fails at the vacuum when limp^op'(p) = 0. The Euler system is also 
genuinely nonlinear in the sense of Lax [38 1 (away from the vacuum), provided the pressure satisfies the 
convexity condition 

p"{p) + pp'{p)>0 (p>0). (1.3) 

These two standard assumptions are made in the present study. (See Lax [38J and Dafermos IT4l for a 
background on nonlinear hyperbolic systems.) 

In continuum physics, compressible fluid flows are often governed by general equations of state p - p{p) 
and, in addition, the selection of shock waves in such flows may be driven by possibly nonlinear, higher- 
order modeling terms. The physical models typically contain second-order derivatives of the unknowns 
and take into account the effect of the viscosity, capillarity, and Hall term in the fluid under consideration. 
It is our objective in the present series of papers to develop the mathematical tools which are required 
for investigating the existence and properties of complex fluid flows. In particular, it has been recognized 
in the past fifteen years that "nonclassical" wave structures and dispersive shock waves do arise in such 
flows, which are not observed with polytropic fluids and with the vanishing viscosity method, whereas 
they do play critical role in many physical applications and, therefore, deserve the attention of applied 
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mathematicians. We refer to LeFloch ||43l for a review of some mathematical tools, and to Section [L4l below, 
for numerous references in a variety of areas: van der Waals fluid dynamics; dispersive shallow water flows; 
quantum hydrodynamic; Bose-Einstein condensates; Boussinesq model; Green-Naghdi models. 

As far as plane-symmetric solutions to the Euler system are concerned, the vanishing viscosity problem 
was first studied by DiPerna in the pioneering work flSl [T9ll . DiPerna relied on Tartar's compensated 
compactness method [58, ,59 J and a compactness embedding theorem by Murat {56, ,57J. This method 
was later extended by Ding, Chen, and Luo |[T7|, Morawetz l55ll . Lions, Perthame, Souganidis Il48l , Lions, 
Perthame, and Tadmor ||49||, Chen and LeFloch | [Tll[T2| . and LeFloch and Shelukhin flS^. In these works, 
all solutions under consideration have arbitrary large but bounded amplitude. 

Several years ago, LeFloch and Westdickenberg 1471 opened the way to constructing solutions within 
the broader class of solutions with finite energy and realized that working within such a large class of 
solutions was necessary in order to overcome certain limitations in DiPerna's theory. In this class, they 
established the first existence result of radially symmetric fluid flows — including the singularity at the 
center of radial coordinates. The present series of papers follows this strategy and, by exhibiting suitable 
algebraic and differential properties of the Euler equations and by covering broad classes of constitutive 
laws, we develop a general tool in order to handle the complex fluid flows arising in physical applications. 
In this first paper, we focus on the zero viscosity-capillarity limit and on flows with plane symmetry, and 
we refer to the companion papers | |24ll25ll for other aspects of the proposed method. 

Throughout, we adopt the following standard notation. We denote by C a positive constant and by 
C{(Sa) a constant depending on some known quantity Sq; the specific value of this constant may change 
from one occurence to another. Finally, it will be convenient to write A < B when the inequality A < CB 
holds for some constant C > 0, while we will use the notation p(p) ^ p'' if p{p) = Cpi'. The so-called 
"Japanese bracket" (x) stands for Vl + (for a real variable x). 

1.2 Navier-Stokes-Korteweg system for viscous-capillary fluid flows 

We consider the Cauchy problem associated with the Navier-Stokes Korteweg NSK system 

pt + (pw),T = 0, 

{pu)t + [pU^ + pip))^ = L[p, U]x + K[p]:c, 

L[p,u]--^i{p)ux, ^"^'^^ 
K[p] := pK(p) P.V.V + 2 (pk'(p) - k(p)) p^, 

with initial conditions 

ip, pu)\t=o = ipo, Pouo) onD (1.5) 
for some prescribed functions po > and Uq defined on D. Here, the nonlinear viscosity = |ti(p) and the 
nonlinear capillarity k - K(p) are prescribed and smooth functions in p > 0, which may be singular at the 
vacuum. Throughout and for definiteness, we always assume that the following (mild) assumptions are 
satisfied: 

(1) pip) > (for all p > 0), 

(2) limu(p) = 0, liminf u(p) > 0, ^'^■^^^ 

p—>0 p— » + oo 

(3) either k = or K(p) > (for all p > 0), (1.6b) 

(4) pp'<p, p\p'\<p, p\k'\<k. (1.6c) 



The (mass density and velocity) dependent variables (p, «) € ]R+ X ]R are unknown functions of {t, x) e 
IR"*" X D. We distinguish between two cases: 

• D = T and the NSK system is thus set on the torus or, equivalently, a bounded interval with periodic 
boundary conditions. 

• D = ]R and the NSK system is set on the real line, subject to the following asymptotic condition, for a 
given constant mass density p* > 0, 

lim (p, pu){t, x) = (p*, 0), t > 0. (1.7) 
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Solutions to jlAl may exhibit cavitation phenomena, in the sense that the mass density p may vanish. 
While this phenomena is unavoidable in solutions to the Euler system jl.li , it can be avoided in jlAi , 
provided the viscosity or capillarity are sufficiently "strong" near the vacuum. While it is straighforward 
to define finite energy weak solutions when the mass density remains boimded away from zero, we will 
need to introduce (in Section|3ll a suitable notion of weak solutions with cavitation in order to handle more 
general viscosity and capillarity coefficients. Observe that if po > we can state the initial condition as 

(p,w)lf=o = (Po/"o) onD. (1.8) 

Consider first the equations l|1.4|l posed on the torus T. The local physical energy associated with the 
NSK system (cf . Section I2l2l below) reads 

E[p, u] := \pu^ + peip) + ^K{p)pi (1.9) 

in which the internal energy e = e{p) is defined by e' (p) := p(p) / p^, while the total energy within T is defined 
as 

^[p,w](0- 1 E[p,u]{t,-)dx. 
Jt 

Interestingly enough, the NSK system is endowed with another energy fimctional, which we refer to here 
as the effective energy. It is based on the effective velocity, defined by 

u:=u + ^p.Y. (1.10) 

More precisely, the local effective energy is E[p, u] - E[p, m] and the total effective energy is defined to be 

I[p,u] ■.= <^[p,u]. 

Observe that both and £ control K^^'^px in L^; and that a combination of S and S controls -^Px in L^- 

Consider next the Cauchy problem dl. 41 - 111. 8b posed on the real line R. The local physical energy needs 
to be renormalized in agreement with l ll.7b as follows: 

1 1 
£*[p,M] pu^ + peip) - P*e{p*) - {e{p*) + P*e'(p*)) (p - P*) + 7:K{p)pl, 

(1.11) 



S-^[p,u\ ■■= Ei,[p,u]dx, 

JR 

while the functionals E* ans tf* are defined similarly. Of course, the expressions given on the torus coincide 
with the ones on the real line by replacing p* by 0. From now on, we therefore use the latter notation in 
both cases. 

Our result about the NSK system concerns the Cauchy problem with initial data {po, Mq) : D ^ [0, +oo) x 
R with finite total (physical and effective) energy, i.e. 

S'Apo, «o] + S'Apo, Wo] < (1-12) 

In the course of writing this paper, we discovered a natural restriction on the capillarity coefficient of 
the Navier-Stokes-Korteweg system, which we refer to as the Strong Coercivity (SC) condition. It turns 
out to play a fimdamental role in the mathematical theory; see Section lZSl below, and especially l|2.26t and 
l|2.27b . Under the condition (SC) and for all finite energy initial data, we will establish the following finite 
total energy-dissipation estimate: 



sup(<#'4p,M](f)+<^*[p,M](o)+ rr ii{p)uidtdx 

t>0 jJ(0,oo)xD 

. [[ i^^^ipxf . '-^p\. . '^^p^\dtdx (1-13) 

JJ(o,c<,)xd\ P P " P I 



/(0,«))xD 
< <f*[po,Mo] + f?*[pO/Mo], 
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for which we refer to the following section. 

Before stating the theorem, we also introduce the no-cavitation (NC) condition 



(NC) For some s > 0, limirtf p'' ^^^|U(p) > or liminf p^'^'*K(p) > 0. (1.14) 

As we show below, the condition (NC) ensures that vacuum does not appear in the fluid if it is not present 
initially. However, the following theorem covers, both, the non-cavitating and the cavitating regimes. 

Theorem 1.1 (Weak solutions to the Navier-Stokes-Korteweg system.). Consider the Cauchy problem ilAi - 
m.SI l associated with the Navier-Stokes Korteiveg system, when the spatial domain D is either the torus or the real 
line and, in the latter case, the condition \\.7\ is assumed for some p* > 0. Assume that the initial data po, wq have 
finite total physical energy and effective energy (see l|1.12|l ), and that the viscosity and capillarity coefficients satisfy 
the strong coercivity condition (SC). Finally, consider one of the following four set-ups: 

(i) The no-cavitation condition (NC) is satisfied and po{x) > Ofor all x e D. 

(ii) D = T; the capillarity coefficient k = vanishes identically; the integral ^pq\uo^^^ dx is finite; for some 
s > 0, the inequality p{p)^ < p^l^p{p) holds for small p > 0. 

(Hi) D = R with p* = 0; K s 0;/or some s > 0, the inequalities p < p(p) < p^ hold for small p > 0; the total mass 
of the fluid ^[po] '■= po dx is finite. 

(iv) D = ]R with p* > 0; the viscosity-capillarity pair (1, k) satisfies (SC); the inequalities p{p) < p^^^ and 
k(p) < ^^-^ hold for small p > 0. 

Then, the Cauchy problem for the Navier-Stokes Korteweg system 1 11 .41 admits a global-in-time weak solution 
{p,u) : [0, +oo) X D — > [0, +oo) X R zvhich has finite total energy and dissipation in the sense ((f^f^). In the 
non-cavitating case (i), one also has p{t,x) > Ofor all {t,x). 

Observe that the right-hand side of the NSK system must be understood in the sense of distributions. 
Theorems l3.1il3^l3.5l and l3.6[ stated and proven below, correspond to the above four cases, respectively. 

While many well-posedness results are available in the literature for augmented versions of the Euler 
system such as the NSK system l ll.4t , as well as variants and multi-dimensional generalizations, most 
papers, however, restrict attention to the purely viscous case (k - 0) and to non-cavitating solutions. We 
do not attempt here to review the vast literature on the Navier-Stokes system and only quote works on 
one-dimensional solutions. Assuming k = 0, Hoff |3ll first treated the case of constant viscosity and finite 
energy solutions. More recently, Mellet and Vasseur 1531 covered the viscosity coefficients p(p) = p" with 
a < 1/2 (a condition which prevents cavitation) and initial data with finite energy and (in our terminology) 
finite effective energy. More recently, for the same class of initial data and viscosity coefficients, Jiu and 
Xin ||33J treated the interval a > 1/2. 

The study of the Euler equations with a capillarity term was tackled much more recently. Under various 
conditions on the viscosity and capillarity coefficients, strong solutions were constructed in Danchin 
and Desjardins IfTSi , Benzoni-Gavage, Danchin, Descombes |3j |4J, and Hao ||27| , while the existence of 
weak solutions was established by Bresch and Desjardins ||6), Jiingel |35l , as well as Gamba, Jiingel, and 
Vasseur |211 . In contrast, the existence theory proposed here relies on the strong coercivity condition (SC), 
only. 



1.3 The zero viscosity-capillarity limit 

Our second objective is the convergence with finite energy solutions to (NSK) in the limit of vanishing 
viscosity and capillarity. We therefore replace p, K in the right-hand side of (11.4b by coefficients p^, 
depending on a parameter e and approaching zero when e ^ Q. For definiteness and without genuine loss 
of generality, we rescale the given viscosity and capillarity coefficients p, k, as follows: 

p^p) = ep(p), k'(p) = 6(e)K(p), (1.15) 
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where b = 6(e) tends to zero with e. Then, we consider the finite energy weak solutions {p^, if) associated 
with the NSK systems : 



p',+{p'u% =0, 
(ph^t + {p'iiff + p(pnl =1^9', 11% + Knp%, 
U[p',u']~p%p')ul, 



(1.16) 



r [p^] -p'K^p') pl, + \ (p'K^'ip') - K'{p')){pl)\ 



and a family of initial data p^, 



{p\p'u%=o = {pl,plul). (1.17) 



We thus consider the singular limit problem e ^ and we are going to establish that finite energy weak 
solutions to the Navier-Stokes-Korteweg system converge to finite energy solutions to the Euler system 

In the case = and restricting attention to polytropic fluids p(p) - p'' (with 1 < y < 3), Chen and 
Perepelitsa flSlI first established a convergence result of the form above; they restricted attention to the 
viscosity coefficient p^{p) = e and, by following LeFloch and Westdickenberg's method |[47). observed 
that the diffusion term can be controled by a priori estimate derived earlier by Kanel 1361 for a different 
purpose. Next, Huang et al. 1321 treated the viscosity functions p'^{p) = ep" and polytropic fluids with 
I < a < y. More recently, Charve and Haspot fTO) were the first to tackle the general viscous-capillary 
problem and established a convergence theorem for polytropic fluids, the viscosity coefficient fi'^(p) ~ ep, 
and the capillarity coefficient K'^(p) = e^p~^. 

In the present paper, we cover a broad and physically realistic class of viscosity, capillarity, and pressure 
functions. Our main restriction imposes that the capillarity term is "tame" with respect to the viscosity. 
In the theory of nonclassical solutions to hyperbolic conservation laws (reviewed in [43]), it has been 
recognized that the capillarity should be bounded by the square of the viscosity, that is, 

K(p) < p{pf, (1.18) 

for, otherwise, oscillating patterns would be generated in the limit of vanishing capillarity and would 
overcome the smoothing effect of the viscosity term. For instance, this inequality can be justified by 
considering the behavior of traveling wave solutions (or by numerically computing the vanishing viscosity- 
capilarity limit; cf. l29l ). When the dispersion effects are dominant, the method adopted in the present 
paper does not apply and the limit, in general, fails to be a weak solution to the associated hyperbolic 
system: this issue was first investigated by Lax and Levermore l39l l40l |4T1 in their work on the small 
dispersion limit of the Korteweg-de Vries equation. 

The inequality l|1.18l l, while being a general feature of vanishing diffusive-dispersive limits of hyperbolic 
systems, is, in the present setup, valid if p is close to a constant, say p ~ 1. In order to take into account 
the scaling of the NSK system near the vacuum, we now consider (for the purpose of motivating our 
tame condition) mass density functions p that remain close to some constant r > 0, say. By rescaling an 
(NSK) solution p by r, we obtain a new solution p = p/r to (NSK), with rescaled viscosity and capillarity 
coefficients 

p(p) = r-^p{rp), K{p) = rK{rp). 

When applied to this rescaled system, the condition l ll.lBI I becomes r K(rp) < ^'^|P , in which p remains close 
to 1, and r > is a parameter. By assuming that the implied constant is independent of r, this inequality 
leads us to the condition 

which is the main restriction required in order to validate the zero viscosity-capillarity limit. The condi- 
tion (|1.191 is thus formulated as 

(TO I ^P^^'-f (1.20) 
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and is refered to as the Tame Capillarity (TC) condition. As explained above, this condition is necessary 
since, when it is violated, highly oscillating patterns arise in solutions to (NSK) and prevent their strong 
convergence as the viscosity and capillarity approach zero (cf. again 1431 ). 

Our second assumption restricts the growth of the viscosity and capillarity coefficients and, specifically, 
we impose the following Growth Rate (GR) condition (for all p > 0): 

(GK) (L21, 

px +5k> 0. 

The second inequality, essentially, requires the lower bound k > p~^, or else k vanishes identically. Recall 
finally that, imder the coercivity condition (SC), the energy inequalities associated with (NSK)f yield the 
following uniform boimd for the solutions p^, u^: 

sup (£ ifm + £4p^ u'm) + f[ ^{pMf dtdx 

f>0 JJ(0,ix))xD 



+ 



(1.22) 



<'^*[pl,ul] + ^^[pl,ul], 



in which we now assume that the total energy of the initial data remains uniformly bounded as e ^ 0. 

Theorem 1.2 (The zero viscosity-capillarity limit). Consider the Navier-Stokes-Korteweg NSK systenie posed on 
the real line M (with p^, > 0) and for polytropic pressure laws p - p'' with y e (1, 5/3]. Assume that and p^ have 
the form il.lSh where the viscosity and capillarity coefficients satisfy the tame condition (TC), the strong coercivity 
condition (SC) and the growth rate condition (GR). Consider global-in-time solutions {p^, u^) : [0, +oo) x R — > 
[0, +oo) X ]R associated with a family of initial data with uniformly bounded total (physical and effective) energy 

\m\swpUi,[pl,ul] + (f*[p^,M^]) < +00 

and satisfying the energy-dissipation inequality il.22h Then, as e — > 0, the family of solutions {p^,u^) converges 
almost everywher^ (after possibly extracting a subsequence) toward a limit (p, m) : R+ x R — > [0, +oo) x R which 
is a weak solution with finite total energy to the Euler system ll.l^ . When p* = and the total mass is uniformly 
bounded at time t = 0, that is, 

limsup ./#[Pq] < +00, 
then the limit has finite total mass for all times, that is, 

t, ■)] < +CO, t > 0. 



For simplicity in the presentation, our convergence result is stated for polytropic fluids and for plane- 
symmetric fluids defined on the real line, although these restrictions can actually be removed and we refer 
to l24ll25l for further details. Finally, from Theorem 1 1.21 we recover an existence result first established in 
ma- 
Corollary 1.3 (Existence theory for the Euler system with finite energy data). Given any initial data (po, Mq) 
with finite total energy 

<^*[po,Mo] < +00 

for some p-i, > 0, the corresponding initial value problem associated with the Eider system of polytropic perfect fluids 
admits a global-in-time solution (p, w) : R+ x R — > [0, +oo) x R, which has finite energy for all times, with 

SAp,um<S^{pQ,UQ\, t>0. 

Furthermore, when p* vanishes and the total mass is initially finite, then 

^[p{t, •)] = ^[po] < t > 0. 



^ At points where the limit p vanishes, the momentum p' u'' converges toward pu but the velocity need not converge. 
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Recall that, if the condition (NC) is assumed, the solutions to (NSK) given by Theorem 11.11 satisfy the 
non-cavitating property infjR p''{t) > (for all t > 0). In constrast, "general" solutions to the Euler system 
given in Corollary 11.31 may always contain vacuum regions. Furthermore, the solutions in Corollary 11.31 
are shown to satisfy only the conservation of mass and momentum, while, in LeFloch and Westdickenberg 
1471 , the weak solutions were proven to satisfy all of entropy inequalities (associated with subquadratic 
test-functions). 

1.4 Shallow water flows, quantum hydrodynamics, and Boussinesq models 

First of all, the Navier-Stokes-Korteweg system l|1.4|l describes the isentropic flow of a compressible fluid 
represented by its mass p > and velocity u e R, subject to a viscous force ^{p)ux, when the internal 
energy admits the decomposition e(p) + K(p)(pj:)^. The constitutive behavior of the fluid is determined by 
prescribing the viscosity function = [^{p) and the capillarity function k - K(p), as well as the pressure law 
p = pip) (or, equivalently, the internal energy e = e{p)). 

However, aside from modeling a compressible fluid subject to viscous and capillary forces, the system 
l|1.4l l also arises in many other physical applications, including in quantum hydrod3mamics and in the 
theory of water waves. Let us present the relevant expressions of the functions p, k, and p for each of these 
model. 

• Polytropic fluids and van der Waals fluids. The pressure fimction is classically assumed to be the one 
of a polytropic perfect fluid p(p) ^ p'>' with y e (1, +oo). In the kinetic derivation of the Navier-Stokes 
system presented in |9], the viscosity coefficient p is a prescribed function of the temperature T, 

specifically - Vt. For a polytropic perfect fluid (with fixed entropy), one thus has T =^ p' which 
leads to the law p{p) ^ pO'-i)/2 for the viscosity. Another classical equation of state which describes 
complex fluids beyond ideal fluids is given by the equation of van der Waals (after a standard 
normalization) p(p) = a p'' (3 - p)~'' - 3 p^ (with p < 3) for some adiabatic exponent y > 1 and with 
a := 8e^'''~^'^^^ (the constant S representing the entropy). 

• Shallow water flows. The Saint Venant model, also called shallow water model, is formally identical 
to the Euler equations jl-l} but corresponds to the pressure law p(p) ^ p^. The viscous shallow water 
equation (for instance derived in |j7| and also in |23ll ) corresponds to the coefficients p{p) — p and 
K = 0. Finally, in order to include surface tension effects as in |52| , we can take K(p) ^ 1. 

• Quantum hydrodynamics and Bose-Einstein condensates. It is well-known that the Madelimg 
transform turns a nonlinear Schrodinger equation of the form idtu - An = f{\u\^)u into the NSK 
system l|1.4|l , in which p - f, p = 0, and K(p) ^ 1/p. A typical expression for the function p(p) is 
the linear fimction p (for the cubic NLS), but a wealth of more refined models exists such as, for 
instance, p(p) = Ap^' + Bp^^; cf. |37| for further details. A closely related model is provided by the 
quantum Navier-Stokes equation, described in Harvey (28l| and fS], which takes the form (|1.4l l with 
the coefficients p{p) - p and K(p) - 1/p. This model is also discussed by Hoefer et al. [30J in 
connection with Bose-Einstein condensates and regarded (after transformation) as an extension of 
Gross-Pitaevskii equation. 

• Boussinesq model and generalizations. It is also instructive to write the system l|1.41 in mass 
Lagrangian coordinates {t, y) defined by dy = pdx — pudt and to introduce the new dependent variable 
17 = 1/p; cf. Section IZTl and equations (|2.7b . As observed in |3|, the Euler equation in Lagrangian 
coordinates coincides with the Boussinesq equation when K(p) = p~^. The standard Boussinesq 
equation corresponds to p(p) = p^, but it is possible to consider more general pressure laws; cf. Bona 
and Sachs [SJ for a class of generalized Boussinesq equations. See also Green and Naghdi [,26] and 
Lannes and Bonneton l46ll . 

The above examples suggest to consider pol5momial (or rational) functions p, p, and k. Although our 
framework is much more general, it is thus interesting to indicate the range of application for our main 
results when p, p, and k are power laws of the form 

pip) = Pop'', pip) = fop", k(p) = Kopl^. 
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First of all, the conditions (NC), (SC), (TC), and (GR) are then equivalent to 



(NC) a < i or /3 < -2, 
(SC) 2a - 4 < jS < 2a - 1, 
(TC) jS = 2a - 3, 
(GR) /3 > -5, 

respectively. Observe that the constraint /j(0) = imposes a > 0. It is assumed throughout Theorem ll.il 
that (SC) holds; furthermore, the four items in this theorem correspond to 

1 

(0 a < - or < -2, 

(ii) Ko = 0; 2y > a, 

{in) kq = 0; < a < 1, 

{iv) - 4 < /3 < -1; > 2a - 3; 

respectively. The assumptions in Theorem ll.2l correspond to 

/3 = 2a - 3; jS > -5; > §' 

respectively, so that the latter two inequalities are equivalent to a > 2/3. 

1.5 The finite energy method for augmented Euler systems 

We complete this introduction with a sketch of the method of proof developped in this paper and the 
companion papers |24ll25ll . Recall that the proposed method is built upon a strategy first introduced by 
LeFloch and Westdickenberg |47J in order to cope with geometrical effects in flows with radial symmetry 
or within a nozzle. The present work allows us to encompass augmented versions of the Euler system for 
real compressible fluids, as well as to study singular limit problems. In the rest of this section, we especially 
emphasize this method for the vanishing viscosity-capillarity problem. 

1. Initial data with finite total (physical and effective) energy for augmented Euler systems. Our 

main assumption is that the initial data of the Euler system jl.li have finite energy, only. When the problem 
is posed on a compact domain, such as the torus T we thus assume that the total energy S'i,[po, Hq] of the 
initial data (po, Uq) is finite. On an unbounded domain such as the real line IR, the asymptotic limit p* > of 
the mass density at infinity must be specified and we must use the normalized energy l|l.lHl . Furthermore, 
when pi = is chosen to vanish and the total mass ^[po] is also finite, then this condition also holds for 
all times. 

Given any augmented version of the Euler system such as the system ||1.4| |, we naturally impose that 
the augmented total energy is finite at the initial time. This energy now takes into accoimt contributions 
associated with the augmented terms and, specifically for the NSK system, the capilarity contributes 
jK{p)pl. Importantly this term has a favorable sign, when the capillarity is positive, as is implied by the 
physical modeling. 

Furthermore, we observe in this paper that the augmented system l ll.4|l also admits an effective total 
energy, obtained by a suitable transformation of the unknowns of the original system. The effective 
velocity ii defined in 111. 41 is introduced and the corresponding effective energy £i,[p^,if]{t) contains the 
term j jj{p)-^ pi / p^ , which is now a contribution of the viscosity and, again, arises with a favorable sign. This 
effective energy was first used by Bresch and Desjardin for the Navier-Stokes system ID. 

2. Global-in-time finite energy solutions to augmented Euler systems. Our first task is to establish 
the existence of global-in-time solutions to the augmented system under consideration, when the initial 
data have the integrability and regularity properties implied by the finite energy condition, which, for the 
NSK system, takes the form l|1.12|l , only. Whenever the cavitation phenomena can be avoided and the 
mass density p remains bounded away from zero, it is straightforward to define a notion with finite energy 
weak solutions in the sense of distributions. Yet, in order to establish the existence of weak solutions, an 
important structure conditions on the augmented terms is required: while the physical energy is naturally 
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dissipative, a condition arises for this property to hold for the effective energy and, in this paper, for the 
NSK model, we introduce the notion of strong coercivity, as will presented and investigated in Section |2l 
below. 

However, cavitation usually occurs in solutions to augmented models when very general constitutive 
laws are considered. For instance, in the NSK model, this is the case when the viscosity and the capillarity 
are too"weak" near the vacuum; see the condition (|1.14|l , above. To handle the cavitation phenomena in 
augmented models, we need a notion of weak solution with cavitation. For the NSK model, this issue is 
discussed in the second part of Section|3l below. 

3. Higher-integrability property of the pressure. Next, in order to analyze the singular limit {e 0, say) 
when the augmented model formally converges to the original Euler system, we assume that the physical 
and effective energies of the augmented model are uniformfy bounded with respect to the parameter e. Our 
first task is to derive several additional higher-integrability properties of the solutions to the augmented 
system, which allow us to get a better control on the solutions. Importantly, the bounds should be uniform 
as e ^ 0. Standard parabolic-type or dispersive-type bound simply blow-up when e ^ 0, and from now 
on we must exhibit additional structure from the Euler system. 

Our first technique extends an argument in De Lellis, Otto, and Westdickenberg IT6ll (for scalar conser- 
vation laws) and LeFloch and Westdickenberg (47| (for the Euler system), and combines the conservation 
laws for the mass and momentum, as follows. For any system of the form 

{p'u')t + ip'iu'f + p{p'))^ = y^[p^ u%, ^^'^^^ 

with augmented terms denoted by X^[p^, u^] and Y^lp^, u^], we can introduce a function by setting 

K := p^ K ■■= -P'"' + Xnp', u% (1.24) 

and then write 

{p'u'h')t + (p'iu'fh' + pip'W)^ = p'p(p') + p'u'X'lp', u'] + h'Y'[p', if]. (1.25) 

This identity is used as follows; cf . Section l4Al below, for the NSK system. By multiplying this identity by a 
positive test-function 6 = d{t, x) and after a suitable integration argument, we obtain the uniform estimate 
on the spacetime integral 

JJ [p'pip') + p'u'X'[p', u'] + ¥Y'[p', u']) e dtdx. (1.26) 

At this juncture, our main observation is that the main term p^p{p^)6 is non-negative, and this approach 
eventually leads us to the spacetime estimate 

pW)eLL (1-27) 

while the contributions from the augmented terms turn out to be controlable by the energy-type estimates 
already established. For instance, for the Navier-Stokes model which contains the viscosity term Y^[p^, u^] = 
p{p^)u^^, we check below that a mild condition on the growth of the viscosity guarantees that the average 
of Y'^[p'^, M*^] is controled by the initial (physical and effective) energy. Furthermore, as observed by LeFloch 
and Westdickenberg [475, this argument applies in radial symmetry and leads to a uniform estimate valid 
even at the center of symmetry. 

4. Higher-integrability property of the velocity. A better integrability property for the fluid velocity 
must be derived next an this is achieved with suitably chosen mathematical entropies of the Euler system, 
following an idea in Lions et al. |49]. Recall that the family of weal|3 entropy pairs is generated by an 
entropy kernel, denoted below hy x - x{PrU;v), and an entropy flux kernel, denoted by = o{p,u;v). 
Specifically, for any continuous function ip - ijj{v) (with subquadratic growth, say) we can introduce 



if 



''{p,u)-- \ x{p,u;v)\l){v)dv, q'^{p,u):= I a{p,u;v)ijj{v) dv (1.28) 



^ That is, entropies vanishing on the vacuum. 
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and, for any smooth solution to an augmented Euler system, derive additional balance laws of the form 

dfT^'^'ip', tf) + ^.,^'''(p^ u') = Z'[p', u']. (1.29) 

Here, the right-hand side Z^[p^, if] vanishes for smooth solutions to the Euler system and in addition, for 
solutions to the augmented model, can also be controled by the energy-type estimates already established. 

The second higher-integrability estimate is now obtained, by integration of l ll.29ll , and by observing 
that we control the entropy flux 

sup r /(p^M^)(^,x)df (1.30) 
X Jo 

in terms of the total entropy and a spacetime contribution in Z^[p^,u^], both latter terms being already 
controled by the existing estimates. Here, our main observation is that, for the family of functions ip - v\v\, 
the inequality (11.291 1 has a non-negative entropy Bux q^'{p, u). The flux in lll.30l l grows typically like \u\^ in 
terms of the velocity. 

5. Young measures for finite energy weak solutions. Equipped with the estimates l|1.27|l and l Il.SOIl we 

are then in a position to introduce a Young measure, say v : 1R+ x ]R ^ Prob(]R+ x R), in order to represent 
all weak limits of expressions like /(p*^, u^). Following [47|, we rely first on the higher-integrability estimate 
(|1.271 for the mass density and establish the weak convergence 

/(p^lO-><v,/>:= rr fip,u)dv (1.31) 



for all continuous function/ = /(p,M) satisfying the growth condition |/(p,M)| < /o(p) pp(p) with limp^+oo/o(p) 
0. Next, we take into account the higher-integrability estimate | |1.30l l for the velocity and we check that we 
can allow a velocity behavior of the form fi{u) \u\^ with lim|u|^+oo /i(w) - 0. 

At this juncture, we point out that the arguments above apply to the Cauchy problem posed on the real 
line, but need some adaptation to apply to the torus. Here, we can also rely on a property of propagation of 
equi-integrability for the velocity, first proposed in {47]. Here, the basic strategy is to integrate the entropy 
balance law (|1.29t with, in l ll.28|l , functions i/' with suitably chosen support in the velocity variable. 

6. Reduction with finite energy Young measures for real compressible fluids. Our next task is to 
derive and analyze Tartar's commutation relation |59] satisfied by the Young measure v for every pair of 
mathematical entropies, that is, for all ipi, 1/^2 we rely on the div-curl lemma and establish that 

df^q'^^ - Z'/]'/'^) = <7]*''0</^) - iq'^W'-) (1-32) 

at almost every point (f, x). The higher-integrability properties above are essential in this derivation, in 
order to allow all fimctions i/'i, i/'2 with subquadratic growth at infinity. 

At this jimcture, a major difficulty is to deduce from l ll.32|l that v reduces to a Dirac mass at each point 
{t,x), at least away from the vacuum, which is equivalent to the strong convergence of p'^ and p'^u^. This 
is done by exhibiting some unbalance of regularity between the two sides of (|1.32t . For polytropic perfect 
fluids, this was done in the references cited above and generalized in Il47l to possibly unbounded Young 
measures with finite energy. The generalization to real fluids is presented in the follow-up paper Il24l . 

7. Global-in-time finite energy solutions to the Euler system. The above steps have thus allowed us 
to fully validate the passage to the limit e — > 0. The relevant notion of a finite energy solutions to the Euler 
equations, first introduced in ||47|, yields that, in particular, such a solution (p, u) satisfies the boimd 

[p(f), u{t)] dx < +CO, t > 0. (1.33) 

In the present work where we take physical viscosity as well as capillarity terms into acount, the entropy 
inequalities need not hold, however. 

8. Subcritical and critical scalings. Our theory covers the regime where the capillarity is dominated 
by the viscosity, in the sense of the tame condition lll.20l l. In the subcritical scaling 

^ ^ 0, (1.34) 

the capillarity terms are "negligible" in the limit and it is expected that the solutions we obtained in 
Corollary II .31 satisfy all entropy inequalities (with sub-quadratic growth in the velocity variable). This 
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property can be checked for traveling wave solutions, at least, along the lines of |[TJ|5T1. On the other hand, 
the most interesting regime from the mathematical and physical standpoints, arises in the critical scaling, 
when the diffusive and dispersive effects within the augmented Euler model are "kept in balance", in the 
sense that 

6(e) = a (a fixed). (1.35) 

Then, dispersive terms generate genuine oscillations which (in the limit e ^ 0) drive the effective dynamics 
of "dispersive shock waves", so that a different selection mechanism may be observed and shock waves 
may fail to satisfy standard entropy conditions llTll2l l34l . However, in the present paper, since the pressure 
function is assumed to satisfy the genuine nonlinearity condition dl.SI l, we again conjecture that the entropy 
inequalities are satisfied by the solutions constructed in Corollary ll.3l Again, this property can be checked 
for traveling wave solutions, at least. 

The rest of this paper is organized as follows. In Section |2l we consider the Navier-Stokes-Korteweg 
system, discuss basic algebraic properties, and introduce our strong coercivity condition. In Section |3l we 
establish an existence theory for the NSK system by constructing weak solutions when the initial data have 
finite energy and the viscosity and capillarity functions satisfy certain mild conditions. Next, in SectionlH 
we establish our two higher-order integrability properties for the Navier-Stokes-Korteweg system and we 
conclude with the strong convergence of weak solutions to the NSK system toward weak solutions to the 
Euler as the viscosity and capillarity tend to zero. 



2 Conservation laws and the strong coercivity condition 
2.1 Derivation in Lagrangian coordinates 

The Navier-Stokes-Korteweg system is derived (in mass Lagrangian coordinates) as follows (cf. f225 for 
details). We denote by (t, y) i-* x(i, y) the so-called mass Lagrangian map, defined so that the integral 

J„ /f (0/ y) represents the total mass which was initially located in the interval [a, b] and, moreover, the 
mass initially located at some point y e D has moved to x{i, y) e ID at the time t > 0. From this map, we 
define the specific volume v = 11 p (or equivalently the density p) together with the velocity v by 

u := xt, V ■■= Xy 

For the sake of simplicity in the notation, we keep the same notation for constitutive functions expressed in 
Lagrangian or in Eulerian coordinates. We proceed by prescribing an internal energy function of the form 
e = e{v, Vy), and we postulate that the following action (on a time interval [0, T]) 

a(x):= rr (e{v,vy)-u^ll)dtdy= ff (e{Xy,Xyy) - X^t l^)dtdy (2.1) 

JJ[0,T]xD JJ[0,r]xD 

is formally extremal among all such maps X- It is easy to derive the Euler-Lagrange equation associated 
this variational problem, namely 

^« + - 4^^Xy,Xyy) + \^{Xy,Xyy)\ ) = 0. (2.2) 



dXy \dXyy 



y'y 



Next, by observing that Vt = Xyt - Uy and introducing the pressure function 

P{V, Vy, Vyy) - -^{V, V y) + —{V, l^y) 



dv ^ \ dVy 

we deduce that the unknown state variables u, v satisfy the following Euler-Korteweg system 

Vt - Uy = 0, 

^ (2.3) 

Ut+P{V,Vy,Vyy)y = 0. 

This system includes the effects of the propagation of waves in the fluid described by the pressure function 
P as well as the effects of the capillarity which is modeled by the internal energy function e. 
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In addition, by prescribing a viscosity function v = v(v), we arrive at the Navier-Stokes-Korteweg 
model in mass Lagrangian coordinates 

Vt - M„ = 0, 

(2.4) 

Ut + P{V,Vy,Vyy)y = {v{v)Uy)y. 

Observe that the local energy E{v, u, Vy) = e{v, Vy) + u^ll satisfies the additional conservation law 

E{V, U, Vy)t + {P{V, Vy, Vyy) w) = (V , V y)) + (v{v) U Uy) " V{v) ul . (2.5) 

\ ... /y \ . ^jj^^ . ,y \ . ,y J 

It remains to comment about the internal energy function. A standard choice made in physics (in phase 
dynamics, in particular) is a quadratic dependency of e with respect to Vy, that is, 

e{v,Vy) = e{v) + A{v)-^, (2.6) 

where A = A{v) is refered to as the capillarity coefficient. Observe that linear terms cannot arise, due to the 
invariance of the physical laws by the transformation y i-^ -y. Consequently, the pressure P splits into a 
function of v and a capillarity term, as follows: 

P{V) = -e' {V), P{V, Vy, Vyy) = p{v) - A' (V) y + {A{V) Vy) y. 

Hence, for the constitutive law (|2.6ll , the Navier-Stokes-Korteweg system takes the form: 

Vt-Uy= 0, 

u, + p{v)y = [v{v) uy)^^, +[A'{v) -f - (A(z;) Vy)j^ 
while the associated energy balance equation reads 

,2 



e(z;) + y + A(z;) y +{p{v)u) 



= [v{v) u Uy) - v{v) My + (w (^y^ vl - (a(j7) Vy)] + Uy A{v) Vy 



(2.8) 



The Lagrangian-Eulerian transformation t h-> y{t,x), defined by yt = -pu and yx = p, allows us to derive 
the Eulerian formulation l|1.4b from the Lagrangian formulation above. It is easy to check that the relation 

vviv) = p{l/v), v^Aiv) ■- Kil/v) (2.9) 

holds between the viscosity and capillarity coefficients in Lagrangian and Eulerian coordinates. 

2.2 Local balance laws 
Conservation law for the mass 

We now record several elementary but fundamental properties of the Navier-Stokes-Korteweg system. 
Observe that the first equation in (|1.4)l . that is 

pt + ipu)x = 0, (2.10) 
simply expresses the local conservation of the mass density. 
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Conservation law for the momentum 

The second equation in l|1.4|l . i.e. 

{pu)t + (pu^ + pip) - p(p)".T - pK(p)p,,. - -(pic'(p) - K{p))pl) = (2.11) 



ut + uu^ + ^-^ p;, = - (fi(p)w.v)^, + (k(p)p.v.v + ^K'ip){p^f)^. (2.12) 



expresses the local conservation of the momentum pn and, in view of the mass equation above, has the 
equivalent form 

P pv— -V.V ■ V---.. ■ 2- 

Observe that all the terms above can be given a conservative form (for instance, for the pressure term by 
introducing the function k{p) := J' p~^p'(p)dp), except the viscosity term which contains an extra factor 
1/p. Hence, in the limit of vanishing viscosity and capillarity and for weaic solution, the expression 
Ut + (m^/2).v + {k{p))x should not be expected to vanish in the sense of distributions. 

Balance law for the energy 

We already introduced the local energy E = E[p,u] of the NSK system by 

1 1 

withe'(p) = Let US now define the local internal forces F = F[p, u] by 

F := -pip) + /.i(p)m.v + pKip)px:, + -(pk'(p) - k(p)) (p^,)^ 

and the local energy dissipation D = D[p, u] as 

D ■- i.t{p)ul. 

Then, the local energy balance law reads 

E, + (uE - uF + K(p)pp.vM.v)^ = -D < 0. (2.13) 
Recall that Dunn and Serrin in Il20ll refer the term K{p)ppxUx, above, as the interstitial work. 

Effective velocity 

Given any w e M, we propose here to define the ly-effective velocity 

u'^' = u^cot^Px. 

By taking a) = 1, we recover an expression introduced first in Bresch, Desjardins and Lin 161 and used by 
Mellet and Vasseur [531, while Jiingel 1351 used a) = 1/2. 

cu-Effective NSK system 

We observe here (cf. the derivation at the end of this section that) the pair (p, u") solves a new system which 
has essentially the same algebraic structure to the one of the NSK system, that is, 

Pf + {pun. = U — Px 



P 

(p u^t + (p(M")' + Pip)) = M"[p, «],. + i<:"[p],, 

M"[p, u] ((1 -cv)p + 0) p,. m"), 

K"[p] := p F"(p)p.o + \ (pF'"(p) - K-(P)) pl 
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which we propose to refer to as the w-effective Navier-Stokes-Korteweg system and in which we have 
introduced the following w-eff active capillarity coefficient 



k" :=K-a;(l-a;)^. (2.15) 

p3 

Observe that the structure of the capillarity terms is exactly preserved, while the viscosity is "split" between 
the mass and momentum equations, and the most important outcome of this transformation is that the 
mass equation has now gained a diffusion term. 

The relevant range for co appears to be the interval [0,1] for, otherwise, the system is not parabolic; 
moreover, when co e (0, 1), this effective system is uniformly parabolic. It is natural also to choose co so 
that k" > — which always holds if co equals or 1, or if &; is sufficiently close to or 1 and the reverse 
inequality in (|1.19l l is assumed. In the present paper, the choice co = 1 will play a central role , and we set 
for the rest of this paper 

M:=iri, k:=1?. (2.16) 

Balance law for the a;-effective energy 

For every a; e [0, 1], we introduce the local effective energy 

E" = E-[p,M] := \p{un^ + pe(p) + ^K"(p)p.?, 
and the local effective energy dissipation 

D" = D"[p, u] := (1 - CO) ^i(p){u^f + cv t^p'{p){p,f + co ^K(p)((p,,)2 + C(p)p^), (2.17) 

with C = C(p) being given by (the rather involved expression below arising first in the calculations) 



3\2 K \p) 



(2.18) 



By a tedious calculation, we can check that 

E';' + E" - u'-'F'' + CO G") = -D", (2.19) 
in which the effective internal forces F" = F'"[p, u] are defined by 

:= -Pip) + il-co) f,(p)M, + pF'(p)p... + ^(pK-'(p) - K"(p)) (p,)2 
and the a)-ef£ective interstitial work by 

G" - ^((pe)' + ^/2)p, + P'(P)PP.VM.V - ^-y-pK{p)p:,P:,:, + ^ ||^J 7c(p) - ^^'(p)) p^. (2.20) 

At this jimcture, we observe that, in | |2.18l l, the function Q need not be non-negative, so that D"[p,u] need 
not be non-negative — in contrast with the physical dissipation D[p, u] which is always non-negative. Our 
"strong coercivity" condition defined below will ensure that D"[p, u] is non-negative in average. Finally, 
when CO is chosen to be unit, we shorten our notation and write 

E:=E\ D:=D\ F :=F\ G:=G\ (2.21) 
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Derivation of an effective NSK system 

The derivation of l|2.14|l from l|1.4) l, in principle, is an elementary matter but may be particularly tedious if a 
brute force calculation is attempted: we propose here a rather natural derivation of this important identity. 
First of all, the treatment of the mass equation is trivial from the definition of the effective velocity, while to 
handle the momentum equation we substract jlAi from l|2.14l l and, therefore, we solely need to establish 
the identity 

(a)p mx)t + {lojpumj + cv^pmfj 

= (conixxp - CO p{u + ojnix)x + co pm'px {u + cvntx)^ - {co{l - co) pn pxx + co{l - co)-{^pn' - n^pfj ■ 

Here, we have introduced the notation ni = m{p) and j^px ='■ nix, as well as m := p^/p^, and we have been 
able to cancel out some capillarity-related terms. By defining a)pm' -: q' and writing 

{a}pmx)t = qxt = qtx = -{q'{pu)x)^ = -{copm'{pu)x)^, 

and then observing that co should be treated as a parameter, we see that the above identity splits into two 
distinct identities 

( - copni'{pii)x)^ + (icopunix)^ = - pUx + pm'px m))^ - (a;( p « pxx + 2 (p ~ ")p?))^' 

{co^pml^^ = (a;^( - ^ m„ + pm'px m.,.))^ + (a)^( p n pxx + i^[pn' - n)p^))^- 

Finally, after removing one derivative in x in each term and getting rid of 00 while splitting the first equation 
into terms that depend or are independent of u, it is sufficient to check the following three identities: 

-q'{pu)x + Ipumx = - p Ux + pm'px u 

= nixxp -pnpxx-^ [p n' - n^pl, 

pml = - P nixx + pm'px mx + pn + - [pn' - n)pl, 
which indeed, in view of our definition of m, n, q, follows easily for arbitrary functions p, u. 

2.3 The strong coercivity condition 

We are now in a position to state several fundamental uniform estimates and state our main assumption 
relating the behavior of the nonlinear viscosity and capillarity coefficients. We are interested in the 
mass and energy equations which do provide us with non-negative functionals (possibly after a certain 
normalization), while the momentum has an indefinite sign and cannot be used to derive an uniform 
estimate. 

The Cauchy problem on the torus 

Integrating the mass equation | |2.10|| in space and time gives us 

■^lp]it)-= I p{t,x)dx= I p{0,x)dx = ^[po]=:.^o, t>0, (2.22) 

which is a constant independent of time. On the other hand, by defining the total energy and total energy 
dissipation by 

<f[p,M](f):= ( E{t,x)dx, 2!{p,u\{t) = ( D{t,x)dx, 

Jt Jt 

we see that the energy equation | |2.13l l yields 

^[p,M](0+ r ^[p,M](s)ds = <f[po,Mo] =:4, (2.23) 
Jo 
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in which D[p, u] and thus !^[p, u] are non-negative. 

Next, by defining the total effective energy and total effective energy dissipation by 



<^[p,u]{t):= I E{t,x)dx, &[p,u]{t)= I D{t,x)dx, 

in which, by jl.lVi , 

D = D[p, u] ■.J^p'{p)pl + t^K{p)(pl, + ap)p^), (2.24) 

we obtain 

<?[p,M](0+ f mp,u]{s)ds = I[po,Uo]=:^o. (2.25) 
Jo 

Observe that, contrary to what happens with the physical energy, the term ^[p, u] in the effective energy 
balance law need not be non-negative. The identity l l2.25ll is going to play a central role in our theory, as 
it provides us with an essential control of the second-order derivative of the mass density, that is, pxx 
— provided we can ensure that the effective dissipation remains uniformly positive. More precisely, in 
this paper we propose the following strong coercivity (SC) condition: for some constant Co > and any 
smooth function p = p{x) > 0, 

(SC) r D[p] dx>Co f I pi, + t^!£M£l dx. (2.26) 
Jt Jt \ P I P 

Clearly, it would be sufficient to assume, for instance that the function C is positive and boimded below by 
p~^ (up to a multiplicative constant), but in fact it is not necessary that C be positive. We refer to Section lZSl 
below, for the derivation of sufficient conditions on p and k guaranteeing that this coercivity inequality 
holds. For instance, it does hold when p{p) - p and K(p) = p" with a e (-2,1). Observe that the term 
^p'{p){px)^ has a different scaling in comparison to the terms and pi, and, although it has a favorable 
sign, we cannot take advantage of it when proving (|2.26l l. 



The Cauchy problem on the real line 

To deal with the problem posed on the real line, we need to introduce a renormalization based on the limit 
p* at infinity and, as already stated in the introduction, we need to introduce 

1 1 
EAp, u] --^pu^ + pe(p) - p*e(p*) - {pe{p))\pi,){p - p*) + -K{p)pl, 

— 1 _ 1 

E*lp,u] :=-p(m)^ + pe(p) - p*e(p*) - {pe(p)y(p^){p - p*) + -K(p)p^. 

Thanks to the hyperbolicity condition | |1.2| |, the function p i-> pe(p) - p*e(p*) - {pe{p)y{pi,){p - p*) is 
non-negative and convex, so that E*[p, u] and £*[p, u] are non-negative. 
Upon defining 

A[p,M](0:= 1 E^[p,u]dx, I^[p,u]{t):= f E^[p,u]dx, 

Jk Jr 

the identities | |2.23| | and | |2.25| | remain valid but the integration domain is changed to IR and the normalized 
energies are used (but the dissipation terms do not need to be renormalized). 

Recall that we use the notation tf* for the torus as well as for the real line, with the convention that 
p* = in the former. The strong coercivity condition (SC) associated with the Cauchy problem posed on the 
real line is stated as follows: for some constant Co > and any smooth function p - p{x) > approaching 
a fixed constant p* > at +oo. 



(SC) 



£D[p]dx>Col[pl.flt^dx. (2.27) 
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2.4 A nonlinear Sobolev inequality 

In the case that and k are power laws, the following theorem provides the key to imderstanding the 
strong coercivity condition proposed in the present work. (Related inequalities can be found in Lions and 
Villani [50].) 

Theorem 2.1 (The strong coercivity condition for power laws). Consider positive functions / : D ^ (0, +oo) 
defined on the torus or the real line and, more specifically: 

• If D = T, then consider f in H^(T). 

• If D = R, then the functions f approach a constant pi, > at ±oo and (/ - p*) e H^(R). 
Then, the inequality 

£ fV^.f dx > i^—^^ f-^if.f dx (2.28) 

holds for any a > 1, in which the constant in the right-hand side is optimal. Furthermore, there is no analogous 
estimate when a = 1, in the sense that if the inequality 

f f{f..fdx>c r f-\f:d^dx (2.29) 
Jd Jd 

holds for all functions f satisfying the above requirements, then the constant c>Oin l|2.291 must vanish. 

Interestingly enough, our coercivity inequality enjoys many symmetries: it is invariant under the 
translation / i— > /(xq + ■)/ well as the multiplication by a constant / i-^ A/ and the dilation / i-^ /(A-) 
(this last property making sense if D = ]R, only). This inequality also depends on our choice of boundary 
conditions, and the conclusion of the above theorem becomes false if, for instance, / is taken to be any 
smooth function defined on the interval [0, 1] (without imposing periodic boundary conditions) or if it is a 
function on R admiting distinct limits at +oo. 

Proof. It suffices to consider smooth functions /, since the general case follows by a straighforward density 
argument. 

1. The torus with exponent a = 1. In order to show that there does not exist c > such that ( |2.29t holds 
for all smooth, positive functions defined on the torus, we argue by contradiction and assume that it does 
hold for some positive c. 

Then, identifying the torus with the interval (-5/ jj, we introduce the family of functions 

Ue{x) = {e + \x\^f^ if|x|<l/4 

(defined for e > and a > 1), which we extend to the torus so that it is smooth, bounded away from zero, 
and enjoys uniform bounds on its first and second derivatives outside of |j. Then, we write 



Q\(UeXdx= aV(e + x2)¥-4dx^ aM'"-'dx = ^i\] , 

J-l/4 J-1/4 J-1/4 Ja-JV4/ 

and, on the other hand, again when e — > 0, 

J (Ue,r) dx -> a\a - If j dx = a\a - 1)- (-J 

If (|229j holds, then 



c I QifQl) dx< I C„,f(Ca,e.vJ dx, 

J-\I2 J-1/2 



and letting e go to zero and using the above gives, for some constant A 



2, ,,2 /IV"-" 20" (IV 
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Finally, letting a go to 1 leads to the desired contradiction. 

2. The torus with general exponent a i= 1. Fix a + -1, and / smooth and positive. Setting h := / A, we 
then observe that the inequality 

J f{f..fdx>c J f-\0dx 

is equivalent to 

J h\h..f dx > ((^)' + c (^)'j h"-\h.f dx. 

Therefore, | |2.28| | follows from Step 1. Finally, the case a = -2 can be obtained by a limiting argument from 
the case a + -2, since the constant does not blow up as fl ^ -2. 

3. The real line. The invariance properties of the equation allow to deduce the results for the real line 
from the results for the torus, at least when is compactly supported. Moreover, an elementary density 
argument then leads to the desired conclusion. □ 

2.5 Application to the strong coercivity condition 

Recall that \i and k are smooth functions mapping (0, +oo) to itself, and define 



:= j^|^p,J (k(p)p,, + ^K'(p)p^) dx. 



(this corresponds to the second term in J D[p] dx, the first one being non-negative if p'(p) > 0). Expanding 
in the above formula and then integrating by parts, one sees easily that 



Jd 



9 XX + ^(P)^ j (2.30) 



with, as was defined in the introduction. 



SU K \(,] p 

Recall that the strong coercivity condition (SC) is satisfied if there exists a constant Co = Co(f(, k) > such 
that, for all function p satisfying the boimdary conditions specified in Section fOl 



(SC) np\>-cJL\^.'^y^^dx. 



The following theorem describes settings where the above inequality holds. Clearly, it is the case when the 
coefficient of (p.v)* is non-negative, but this may not be a physically realistic condition, so we consider the 
following broad classes of functions. 

Theorem 2.2 (Sufficient conditions ensuring the strong coercivity condition). (/) Yix [i, k (positive, smooth 
functions on (0, +oo)), and po > 0. Then, there exists e > such that (SC) holds as soon as \p{x) - po| < e. 

(ii) Assume j.L{p) = p" and K(p) = p^\ Then, the condition (SC) holds for some Cq if and only if 

2a - 4 < /S < 2a - 1. 

If this condition is not satisfied, the functional !>[■] is not even positive. 

(Hi) Assume ^{p) = p" and K(p) = p'^ and that the above inequality is not satisfied. Then, there exists c e (0, 1) 
such that the condition (SC) holds for all functions p satisfying \p{x) - po\ < cpo (for all x and some po > 0). 

(iv) If the function C ^'s positive and there exists cq > such that Q{p) > CQp~^, then the condition (SC) holds. 
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Proof, (i) Set p = po+ £<p, where \(p{x)\ < 1 for all x. The functional D becomes 

Taking e sufficiently small, n{p)K{p)/p becomes bounded from below, and the result follows then if D = T 
by the Sobolev embedding theorem, and if D = ]R by the Gagliardo-Nirenberg inequality 

\\<Px\\l < ll</'l|oo||(p.r.vll2- 

(ii) For - p" and K(p) - p^', the expression I l2.30l l becomes 

= ^ [p"^^-^p..f + - 1) + - - 2)) p"^''-Hp.f) dx. 

Comparing this with (|2.28l l gives the condition 

-|S(/5 - 1) + -(a - l)(a - 2) < | ^ 

which results in the range for j8 in the theorem. 

(iii) This follows from the item (i) and a scaling argument. 

(iv) This is immediate. □ 

We conclude with a proposition which suggests a more general characterization of pairs {fi, k) ensuring 
(SC). 

Proposition 2.3 (A characterization of the strong coercivity condition). When p and k are power functions, 
say p{p) = p" and K{p) = p'^ the following three statements are equivalent: 

(i) D[p] > Ofor any positive and smooth function p. 

(ii) 2a- 4 < fi <2a - 1. 

(iii) 1>[p] can be written in the form 



in which F is a smooth and positive function. 

Proof. The equivalence of (i) and (ii) was already established in Theorem l2.2l It is clear that (iii) implies (i). 
So we are left with proving that (ii) implies (iii). Comparing the expression in (iii) with (|2.30l l, one sees that 
the functions 

pip)K{p) F'(p) 
''^P^=^^T(pj' P^'' 
has to solve the ordinary differential equation in p 



If p and K are powers of p, this equation becomes 

ju'(p) - ^p'-"-Pu^ = (i(a - 2)(a - 1) + 1|S(|S - l))p"^P-^ 
Try the ansatz u = Ap^ . First, one needs A = a + /3 - 2. Next, A has to solve the quadratic equation 

- (a + iS - 2)A + (a - 2)(a - 1) + \m - 1) = 0. 



20 



This equation has real solutions if and only if its discriminant is non-negative, that is, 

A = (a + - 2)2 - 3((a - 2){a - 1) + ijS(jS - 1)) > 0. 

Solving this inequality gives the condition (ii). If it is satisfied, we can come back to F, and obtain the 
solutions 

F = p^ s=^(fl + /3-2). 



3 Finite energy solutions to the Navier-Stokes-Korteweg system 

3.1 Existence theory for non-cavitating solutions (k = or > and D = T or R) 

We begin by establishing an existence theory under the non-cavitating condition (NC) introduced in l(1.14|l . 

Theorem 3.1 (Non-cavitating finite energy solutions to the Navier-Stokes Korteweg system). Consider the 
initial value problem associated to the NSK system l|1.4|l posed on a domain D, ivhich is the torus or the real line (with 
p* > fixed), and assume that the non-cavitating condition (NC) and the strong coercivity condition (SC) hold. 
Then, given any initial data {p^, uq) with finite total energy and ejfective energy Sq, satisfying the non-cavitating 
condition po{x) > Ofor all x e D, there exists a global-in-time weak solution (p, u) : [0, +oo) x D ^ (0, +oo) x R 
satisfying the energy estimate {Si,^) with, furthermore, p{t, x) > Ofor all (t, x), and 

floe xloc floe xloc 

where H^, if - denote the standard Sobolev spaces. 

Proof. 1. Reduction to the pre-compactness of solutions to (NSK). We need to establish that the pre- 
compactness of the family of NSK solutions satisfying (at the initial time) a uniform energy bound and a 
lower bound on the mass density. In other words, from any sequence satisfying these uniform bounds, we 
can extract a subsequence that converges (in a suitable topology) to a solution to ||1.4|| . This is the property 
required in order to deduce the strong convergence of approximate solutions to l|1.4|l and, eventually, 
establish the existence of actual solutions. On the other hand, several standard methods are available for 
the construction of approximate solutions; one can for instance use Galerkin-type schemes and we refer to 
Gamba, Jiingel, and Vasseur ||2T1| for further details. 

2. Energy bounds. Motivated by the observation above, we now consider solutions (p", m") associated 
with some initial data (p^, «g) whose physical energy S"^' and effective energy (o^ are uniformly bounded, 
that is, 

sup(^o« + 4") 



< +CO. 



In view of the previous section and our assumption (SC), in particular, this implies the uniform boxmd 
(<f;^)for (p",M"). 

3. Lower bound on the mass density. Here and in the following step, we rely on arguments first used in 
Hoff L31j and Mellet and Vasseur (53] in order to cope with the viscosity term in Navier-Stokes equations. 
Starting with the case of the torus D = T, we observe that the energy provides us with the boimd 

K{p"Mf e LTLi 
while the energy and the effective energy (taken together) control 



(F(P"))^ 

{p"f 



ip'if e L-Ll 
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Thanks to the no-cavitating condition (NC) (stated in (|1.14l l), the latter two estimates imply the existence 
of a > such that {{p")~")x e L^L^. On the other hand, the mass := Jj p" dx is conserved in time and is 
uniformy bounded thanks to our vmiform energy boimd, so that the following lower bound holds: 

supp{t,x)>—. 

.veT I ^ I 

In turn, by Sobolev's embedding theorem, we arrive at the lower bound inff ^ jo" ^ pmin 

> for all n. 

In the case of the real line D = IR with > 0, the same arguments yield again ((p")~").ir ^ ^T^x 
then observe that the bound on £" implies the existence of constants R,r > such that, for any Xq, there exists 
X e [xq - R, Xq + R] such that p"{t, x) > r. We can then conclude again by Sobolev's embedding theorem. 

In the case of the real line D = R with p^, = 0, the same approach can once again be followed. First 
one shows, for any t, the existence of reals R{t) and r(f), depending only on the data, such that p(xo, t) > r 
for some |xol < R- Then one applies the fact that {{p")~")x e i-^i-x to get a lower bound (which goes to zero 
as X — > CO. We skip the details. One obtains, for any T, the existence of a positive function Cr such that 
p"{x,t) > Cr(x,Ofor f < T. 

4. Upper bound on the mass density. In the case of the torus D = T, observe as above that the energy 
and the effective energy (taken together) control ^^'l^p„yi (p")^ ^ ^T^x- Thanks to the mild assumption (|1.6|l , 

this implies that xip") a/"^ (where x cuts off smoothly to p" > 1) belongs to L'^W^' . Indeed, by the 
Cauchy-Schwarz inequality and the growth condition on the pressure, we have 



r ^(p")5.vV^(p^dx< r xip" 

Jt Jt 



Px 



dx 

,,n2 \l/2 



(p«)3/2 - 

J^Xip")e{p")dx J^Xipl^^dx"^ < C(4,4). 



On the other hand, the mass equation gives 

^" 

infp"(f,x) < — ^. 

.V6t'^ ^ ' \T\ 

These two observations imply the existence of the uniform upper bound supj ^ p" < pmax uniformly in n. 

In the case of the real line D = ]R, the second observation above needs to be modified, as follows: by 
virtue of the boundedness of the renormalized energy, there exists a constant pmax > such that, for any 
xq, there exists x e [xq- 1, xq + 1] such that p"{t, x) < pmax- 

5. Uniform estimates in Sobolev spaces. The uniform estimate (S'^^) and the upper and lower boimds 
on p" imply (with uniform bounds) 

II T oo rrl pi 7" 2 r_j2 71" <= f °° f 2 n f ^ F-f^ 

P ^ ^Hoc^xloc ' ' ^floc^-vloc' " ^ '"floc^iloc ' ' ^floc^xloc- 

Furthermore, by returning to the system (NSK) satisfied by (p",u"), we also deduce uniform bounds for 
the time derivatives 

P'l^^nocH-xL <e4oc<oc- 

6. Passage to the limit. Using standard Sobolev compactness theorems and, in particular, Aubin-Lions 
lemma, we deduce that there exists a limit (p, u) such that, for all e e (0, 1), 

weakly-* in L- ^^^^, 

weakly in L^Hj. and strongly in LfH^"*^ locally for all p e [1, +oo), 
weakly in L^H? and strongly in l^ti^'^ locally if k > 0, 

„ J weakly-* in L^L^ locally for all p < +oo, 

I weakly in LfH!^ locally, and strongly in LfH^~^ locally. 

It is immediate to now pass to the limit in the weak formulation of the system (NSK) and check that (p, u) 
solves the equations in a weak sense. 



22 



Observe that, in the case that the capillarity k vanishes identically, then the convergence property 
p" ^ p (weakly in L^H^) is "lost" but, simultaneously, the third-order terms in our system are gone, so that 
the convergence properties above are strong enough to allow us to conclude. □ 

3.2 Existence theory for cavitating solutions (k = and D = T) 

We are now interested in cavitating solutions, which we will obtain as limits of the non-cavitating solutions 
in Theorem l3.1l We begin by considering the Navier-Stokes (NS) system 

pt + {pii)x = 0, 

, , X (3.1) 

{pu)t + [pu"- + p(p)j^ = (pip) M.v).v, 

and we focus on the Cauchy problem posed on the torus (while the following two subsections concern the 
real line). 

Theorem 3.2 (Cavitating finite energy solutions to the Navier-Stokes system). Consider the system i3.1i posed 
on the torus with prescribed initial data {pQ, Uq) with finite energy Sq and finite effective energy Sq. Assume that the 
data satisfy 

polMor^eLi(T) (3.2) 

for some s > and that 

< 1 for small p. (3.3) 

Then, the Cauchy problem associated with admits a global-in-time weak solution (p, u) : [0, +oo) x T — > 
[0, +oo) X ]R, which may contains vacuum regions (and is understood in the sense ofDefinition \3.3\ below), and has 
finite energy and dissipation in the sense {<§-i,2!) with, furthermore, p e L^^ and p{p) e ej(*. 

We emphasize that the mass density p is continuous, so that the following notion of weak solution 
applies. 

Definition 3.3. Under the assumptions and regidarity properties in Theorem \32[ a weak solution with cavitation 

to the system jS.li is defined as follows: the uniform estimate ((f*i^) controls e{p) in L^Lj, pu^ in L'^L\, and p{p)u\ 
in L-^L^, so that the conservative variable and flux terms p, pu, pu^, p{p), and p{p)u\ are locally integrable functions. 
For instance, for p{p)ux and by Cauchy-Schwarz, p(p)ux ^ p{p)i'l + l^ip)- 

Before we proceed with the proof of Theorem 13.21 we establish a preliminary lemma, which takes 
advantage of the viscosity term in order to derive a better integrability property for the velocity. Observe that 
the alternative approach in Theorem l3.61 below, relies on a hyperbolic-type property of higher-integrability, 
which eventually allows us to remove the additional integrability condition (|3.2b . 

Lemma 3.4 (Higher-order integrability of the velocity). Under the conditions in Theorem 1X21 for each T > 
there exists a constant C(T, Sq, Sq) > such that any non-cavitating solution of the Navier-Stokes system \3.1) 
defined on the torus T satisfies for allt <T 

r {p\u\^^'){t,-)dx+ r r \uYp{p)\ux\^dtdx< r po\uo\^^'dx+c{T,so,so). 

Jt Jo Jt Jt 

Proof. We follow an argument in Mellet and Vasseur Il53l and multiply the momentum equation in (|3.11 by 
\ufu. After integrating in space, we obtain 

I {put + puux)\u\hidx + I p{p)x\u\^udx - I {p{p)Ux)x\u\hidx = Q. 
Jt Jt Jt 

The first and third terms above can be rewritten in the form 

p\u\^^'dx + {s + \) f p{p)ul\ufdx, (3.4) 
A + sat Jy Jt 
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while, by integrating by parts and using Cauchy-Schwarz inequality, the second term can be bounded as 
follows: ^ 

r p{p)x\ufudx < f p{p)\uf\u-c\dx < e f n(p)ul\uf dx + - f ^^^\ \uf dx. 
Jt Jt Jt ' ^ Jt V-Kp) 

The first term in the right-hand side above is controlled by the dissipation term in | |3.4| |, provided e is 

sufficiently small. For the second term, we write the integrand in the form p'^^'^if ^Q,y,si2 and we observe 

that: 

(i) The factor p'^^^u'^ is bounded in L'^I^J^ (in terms of the physical energy) and thus in L'^L\. 

(ii) The factor Ji^p^psn belongs to since ll3.3l l is assumed and the mass density is uniformly boimded in 
L^^ (cf. the proof of Theorem lS.ll or see I l3.5b , below). 



Proof of Theorem I3l2l 1. Approximation from the non-cavitating regime. Consider a viscosity coefficient 
- ii{p) and initial data {po, Uq) be as in the statement of the theorem, and let us solve the NS system with 
non-cavitating data (po + 1/^/ Mq) and viscosity /./"(p) = /i(p) + 1/n. By Theorem 13. 11 there exists a solution 
(p", m") which enjoys uniform physical energy and effective energy bounds and such that 



inf p"{t, X) > p'^^ > 0, sup p" {t, X) < praax, 



(3.5) 



in which pj^^ > may depend on n but pmax is independent of n. We are going to establish that these 
solutions converge (as « ^ oo) to a solution to the NS system, denoted by (p, u). 

2. Uniform bounds on the mass density and its derivatives. Thanks to the energy bound, p„ is boimded 

in L'^V^ locally, so that, upon selecting a subsequence if necessary, we assume (n — > +oo) 



On one hand, the inequality 



p" -> p weakly- * in V^L[. 



implies that llfi"(p").i:llL~L5 ~ C((?0/ <^b/ ^o), whereas, on the other hand, 

i(^,"(p"))fi < F"(p")i"vi + ppiip?i Vp^|""|' 

implies that ||p"(p")t|li,2L2+LrL] ^ C(4,4,^o)- 

From the above estimates, we now deduce that f/"(p") is boxmded in 6^^*, indeed: 



|^,"(p")(f,x) - f<"(p")(s,y)| < |f"(p")(f,^) - f'"(p")(f,y)| 



f"(p")(f,y)-:^ r f'"(p"(f,- 



)) 



^" J[y-fc;y+;i] ^" J[y-h;y+h] 

inp"){s,y)-^ r 

^" J[y-h;y+h] 



[y-h;y+h\ 
)) 



^(P"(s,-)) 



upon choosing h = ^J\t - s\. Finally, thanks to Arzela-Ascoli theorem and after taking a subsequence if 
necessary, we have the strong convergence property 



^(p")^^(p) ine|„^^^ 
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3. Uniform bounds on the pseudo-momentum wzji' := ju'^{p")u" , with j.i'^{p) := min(/.i" (p),p). Since the mass 
density is uniformly bounded, it is clear that 

\\m'l\k^L2 < 4, (3.6) 
while its x-derivative can be estimated as follows: 

< ii;{p")K\ + \p:\\u"\< ^"(p")K\ + \p'i\^/p-n \u\ (3.7) 

implying IKm^vllL^Ll. - C((?0/ <^0/ ■'^o)- To estimate the time derivative, we use the NS equations to replace 
the time derivatives of p", u", as follows: 

(mp,=^J(PnK + (f'P'(p"K"" 

= - ii'iip'yi'hii - J_(p(p«)), + J_(^,"(p«)<), - {^-y{p'')p"u'x - {ii;np")p>"u'' 



i.i"{p") 

= ( - (r(p"))x + {^—^i"{p"Kh - Ui;(p"){u"fh) 



+ ( - {^y{p")p'i^"{p"K + ip;yip")p'Xu" - p';(p")u"u^ =■. i + n, 

where the function T = T{p) satisfies r(0) = and T' - t^^£^J£}._ it is now easy to see that ||I||^2 ^-i < 
C{S'o, (S'o,M) while ||II||j^2 < €{£0, S'o, ./#o)- Therefore, by Sobolev embedding, we find 

t loc ^' 

<C(4,4,^o). (3.8) 

H Hoc 

4. Convergence of the fluid velocity. The Aubtn-Lions lemma, combined with the estimates (I3.6|l , I I3.7II 
and (13.8b , gives the existence of that a subsequence of m'^ converges in L^L^ locally to m^. Upon taking a 
further subsequence, we can assume that it also converges almost everywhere. We then set 

[0, p = Q. 

With this definition, w" converges to u almost everjrwhere on the set |p > o|. Next, choose e > Q. Recall 
that /.i(p") converges to ^{p) in 6^^*. In particular, for all sufficiently large n, we have fi(p") > § on the 
set Oe = [pip) > e}/ which gives us a uniform bound for u" in Lf^{Oe)- Taking a further subsequence if 
necessary, we can assume that u" converges weakly to some v in L^^(Oe). It is easy to check that v = Ux. By 
a diagonal argument, we can further achieve that, for each e, m" converges weakly to iix in ^(Of). 

5. Passage to the limit in the momentum. Denoting by x the characteristic function of the interval [-1, 1] 
and for arbitrary L, T > 0, we have 



Ip"m" - pu\ dtdx ^ J \p"^"^{t') ~ ^^'^^ 



The first term in the above right-hand side converges to as n 00 (by the dominated convergence 
theorem), since we have imiform bounds on the integrands. For the second and third terms, we rely on 
the uniform bound on the physical energy and write 



limsup I I \p"u" - pu\dtdx 
)i^+oo Jo Jt 



< c(4)r 
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Letting L tend to infinity, we see that p"u" converges to pu in Lj ^. 

6. Passage to the limit in the momentum flux. The momentum flux is treated in a similar way to the 
previous step, by now writing 

£ \p"{u"f - P"'\ £ J |p"(u")'a- (y) - po'A- ( j)| ltdx 

*l7/"<""'1'-^(¥))'"''"X7/"1'-^(i))''*- 

At this jimcture, we rely on the higher-rntegrability property established in Lemma |3^ and obtain 



limsup I I \p"{u")^ - pu^\ dtdx : 

n^+ca Jo Jd 



C(4)T 



By letting L tend to infinity, we conclude that p"(M")^ converges to pu^ in Lj ^, locally. 

7. Passage to the limit in the viscous force. This final term is more delicate if one wants to cover the 
cavitating regime. We will prove that ii"{p")u" — > ^{p)Ux in the sense of distributions and that all the terms 
involved can be defined. Pick up any smooth function x : ^ [0,1] that vanishes identically on the 
interval [-1, 1], and is identically equal to 1 on [-2, lY- Then, if is any test-function and e > Q, we have 



^ cp^"{p")uUtdx ^ JJ cj,^"{p")u'ixipp-^ dtdx + JJ (/)p"(p«k|i-, 



dtdx 



= I + IL 



Since ju"{p") converges to j.i{p) in G^^'^; for all sufficiently large n, the functions x(^~T^) supported on 

Oe/2- Using that u" converges to in Lf^{Oe) whereas x{^-j-^) H-"{p") converges to x{^) H-"{p") in 6^^*/ 
we have 



dtdx 



when n —> +oo. 



On the other hand, we have 



|n|< 



P"ip") 



Ll. [0,2f] 

Combining the estimates for I and II yields 

\dtdx 



< sup V^ll Vf'"(p")w.?|L2 < sup 

'■' I0,2e] 



JJ (pfi"ip")u",dtdx^ JJ (p^{p)u^dtdx. 



8. Final conclusion We have established that, as m — > +oo and in the sense of distributions, p" — > p, 
p^u" — > pu, p"{u")^ — > pM^ and p"{p")u" — > n{p)Ux. Furthermore, by using the upper bound on p", it also 
follows that pip") — > pip)- Therefore, we can pass to the limit in all the terms involved in the Navier-Stokes 
equations and we conclude that the limit (p, u) is indeed a weak solution. □ 



3.3 Existence theory for cavitating solutions (k = and D = ]R) 

We now generalize our analysis to the real line. 

Theorem 3.5 (Cavitating finite energy solutions to the Navier-Stokes system). Consider the Navier-Stokes 
system i3.1^ posed on R with prescribed initial data (po, mq) with finite energy Sq, finite effective energy Sq, and finite 
mass .Mo- Assume that, for some a > 0, 

p" > pip) > p for p small. 

Then, there exists a global-in-time solution (p, u) : [0, +oo) x T ^ [0, +oo) x ]R fo l|3.Hl , which possibly contains 
vacuum regions and such that {S-k&) holds with, furthermore, p e U^^ and pip) e C]^^. 
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The weak solutions with cavitation above are understood in the sense of distributions, along the lines 
of Definition |3.31 above. 

Proof. Like Theorem |3.21 this theorem is established by approximation from the non-cavitating case treated 
in Theorem 13.51 We will not repeat herre all the steps already described in the proof of Theorem |3.2[ and 
we only emphasize here the novel argument required in the present proof, that is, the higher integrability 
property for the velocity. The following argument relies on a notation which will be introduced at the 
beginning of Section l4!2l below. 

1. A new set of entropies. Given a parameter a e [0, 1], we consider the two entropy-entropy flux pairs 
obtained by choosing the ftmction ip to be [vf'^^ and v\v\" in l|4.7t , that is, 

^-:=,fr\ ^«:=qW", ^:=r,*l", ^:=^*l". 

Using | |4.5t - (|4.71 , these entropies can be checked to satisfy the following pointwise bounds: 

< p^-'+y^+i + p\ur\ \if\ < p(«+i)0+i + p\ur\ 

|rt|<p«^ + |Mr, |7f„| < p"^ + N«, (3.9) 



In particular, we observe that q" is positive. 

2. A bound via the entropy dissipation of rj". We will prove that imder the assumptions of the theorem, a 
solution of Navier-Stokes satisfies 

u 



r r p^"~^^''{^)''~\i{p)u\dtdx<C{T,J^o,So,. 
Jo J p 



Start by writing the conservation law for the entropy rj": 

Wt + {q"h = nMp)u.).. 

Integrating it in time on [0, T] and in space over R gives 

{p{p)Ux)x dx, ds 

^ ~ Jo I '^'""^'^^^"^ ^^'^^ ~ X / '^'"p'^^^'^P^"^ 
where the second equality follows by integration by parts. Using the bounds | |3.9|I leads to 



J p^'-^^\^Y-^ix{p)uidtdx< J (p™+i + piuri)(r)rfx+ J (p™+i + piuri)(o)dx 



-'0 

+ 



The two terms in the above right-hand side can be bounded by using the observation that that p("+i)^+i + 
PImI""*"^ ^ P + pu^ + p''; for the third one, we use Cauchy-Schwarz inequality. This gives 

^ J p^''-^^\^)''-^^{p)uldxdt<SQ + J(Q + ^^ J ^(p)uldtdx + J ^p^^pldxdt. 



lo J P" ■ Jo J ' Jo J P" 

The last term in the above right-hand side remains to be controlled. Since /./(p) > p, it is controlled by 

p3 



/ / ^(P-^-)^ 'i^'^t- therefore get 



r r p^''-'^''{^r-h<p)uUtdx < 4 + ^0 + 4<T), 
Jo J P 
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which is the desired result. 

3. A bound via the entropy flux of rf . We prove here the uniform bound 



r r ipM^^" + p''"""^) dtdx < C{K, Mo, <^o, <^o, T). 
Jo Jk 



which is the desired higher integrability for the velocity. The proof is actually parallel to that of Proposi- 
tion l4.31 the difference being that, here, e = 0, which greatly simplifies the estimates. Thus we can skip the 
details and refer to the proof of Proposition l4.3l Begin by writing the conservation law for the entropy rf 

dtjfip, u) + d^cfip, u) = rf„,{p, w)(fi^(p)M.O.r 

I n HI 

Multiply this equation by sgn(x - y) and integrate it on (t, x, y) e [0, T] X ]R X K. 

By H3.9II , the contribution of I can be boiinded by ■ ■ ■ < C(<fo/ ^o)- Still by (|3.9t , the contribution of II is 



We are left with III. Its contribution can be bounded with the help of the bounds ll3l 
Schwarz inequality: 

I I I sgn(x - y) III dy dx ds 

Jo Jr Jk 

< f^(^K'^P'" + + f'(p)P'""Vxll"vl + ^lip)p^''-'^'{^r'ulj dtdx 



and Cauchy- 



Jo Jr 



2a0 



III, 



III, 



III, 



III4 



III, 



dtdx 



We examine each of the terms appearing above: IIIi is obviously bounded by the energy; III2 was controlled 
earlier; III3 was already treated in Step 2; and the term III5 can be controlled as in Proposition |4.3l All in 
all, 

J (IIIi + III2 + III3 + Illg) dx dt < C{<^o, Sq, T, K). 
Finally, Holder's inequality gives 

^ J IIl4dxdt<|j^ J fi{p)^'''\u\^dxdt^ ~(X / Pl"!^^^'^^) -'^'^0, 
where the prior-to-last inequality follows from the assumption ji{p) < p". □ 

3.4 Existence theory for cavitating solutions (k > and D = R) 

We now consider the fuU system with both viscosity and capillarity terms. 

Theorem 3.6 (Cavitating finite energy solutions to the Navier-Stokes-Korteweg system). Assume the viscosity 
and capilarity coefficients ji,K > satisfy the strong coercivity condition (SO and, that additionally, the pair (1, k) 
also satisfies (SC) and the following growth conditions hold: 



K(p) < 



for small p > 0. 



(3.10) 



Consider the Navier-Stokes-Korteweg system \\A\ posed on the real line with p* > with prescribed initial data 
(po. Mo) with finite physical energy Sq and effective energy Sq. Then, there exists a global solution of(NSK) (in the 



sense of distributions) such that {Si,2i) holds with, furthermore, p e U^^ and p.{p) e nej(*. 
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Definition 3.7. Solutions of(NSK) are understood in the sense of distributions. As already explained in the case of 
Navier-Stokes, the uniform estimate (A^^) ensures that p, pu, pu^, p{p) and ii{p)Ux are locally integrable. It thus 
simply remains to give a meaning to the capillarity term K. It suffices to write it 




and to notice that, since p is bounded, the finiteness of S implies that pK{p)px and (pK'(p) + 3K(p)j pi are locally 
integrable. 

Proof. 1. Approximation by the non-cavitating case. We replace the initial data by (po + 1/n, uq) and we 
replace the viscosity coefficient by [i{p) + l/n. This gives a sequence of solutions (p", u") which, we would 
like to show, is compact and converges to a solution of the desired equation. It can be done in a very similar 
way to the proof of Theorem 13.61 we only describe here two new technical ingredients: the derivation 
of improved integrability for the velocity u (which was proved in Lemma 13.41 in the purely viscous case, 
provided additional integrability already held initially), and the passage to the limit in the capillarity term. 

2. Higher-integrability property for the velocity. We only treat here the case p* > 0, the case p* = can 

be dealt with in a similar fashion. We wish to prove that, for any compact K, J p"|u"p dtdx is uniformly 
bounded in n. This is quite close to the statement of Proposition l4.3l (with e - 1), and indeed the proofs are 
almost identical; note that the assumption (|3.10l l is needed here. The only difference is a new term in the 
last line of (14. 9t , below, specifically 



We indicate here how to treat this error term (following an argument in ||T31 ) and we refer to Proposition l4.3l 
for the rest of the proof and the more general context. Observe first that the bound on S" implies the 
existence of RiSo), r{S'o) such that: for any t, xq, there exists a set I of length one, contained in [xq -R,xo + R] 
such that p{t,x) > r on I. Still by the bound on this implies that there exists R{(S'o),L{S'o) such that: for 
any t,xo, there exists x e [xq - R,xo + R] such that |M(f,x)| < L. Now, pick any xq and find an associated x 
with this property. Then 



3. Passing to the limit in the capillarity term. Recall that, as was established in Theorem l3.2i p" converges 
to p in ej(* away from {p - 0}. Taking advantage of the bound on pl^ which becomes available since k > 0, 
it is easy to show, via the Aubin-Lions lemma, that p" converges to p in L^Hj away from {p = 0}. 

The capillarity term reads K[p]x, with K[p] = pK{p)pxx + \ {pK'{p) - K(p)) pi. We want to show that K[p"] 
converges as n ^ oo to K[p] in the sense of distributions. We focus on the hardest term and show that 
P"k{p")p"x converges to pK{p)pxx ■ Write first 



Once again, we focus on the most difficult term, namely (px)' {p'^){px)^ . Taking x to be the indicator fimction 
of the unit ball, and K a compact set. 





This U° bound on the velocity enables us to treat the term we are discussing: 




p''K{p")p:x = ip"Kip")p:)x - ipKY{p"){p:f. 




< 




Jo Jk \ 



r r x[^)<p)pidtdx 
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The first term in the above right-hand side goes to zero as n goes to infinity, since p" converges to p in L^H^ 
away from {p = 0}. As for the second term, it can be bounded using Holder's inequahty and | |3.10| | 



pe[0,e] 



< sup p'-^ 4), 



in which the upper bound now tends to zero with e ^ 0, thanks to the assumptions Jc(p) < ^^^p- and 
/i(p) ^ as p ^ 0. Finally, the third term can be bounded similarly since for n big enough, p" < | if p < e. 
Therefore, as w +oo 



J /I 

Jo Jk 



Ik 

which is the desired result. 



K(p")(pD' - <p)pi\ dtdx ^ 0, 



4 Finite energy solutions to the Euler system 
4.1 Higher integrability property for the pressure 

We consider in this section solutions of the system (NSK)^ satisfying the uniform bound l ll.22ll , where we 
recall that 

lu'ip) = efi(p), 7c^(p) = 6(e)K(p). 

For simplicity in the notation, we drop the superscript e and simply write (p, u) instead of (p*^, u*^). In this 
section, we assume the following: 

• The equation is set on ]R. 

• The conditions (TC) (defined in (LM ) and (SC) (defined in lIZlTl l) hold. 

(v-l)^ 

• The fluid is perfect and polytropic, that is, p(p) - ^ ' p'' with y > 1. 



As it is usual, we set Q :- '-j-. 



r-i 

2 ■ 

Lemma 4.1. Under the above assumptions, there exists a constanf^ C{(oo, Sq) > such that 

eMpKp)'^C(^o,4) (4.1) 

zvith 

^^1^-^(4,4)0(1), (4.2) 
ivhere o(l) denotes here a function tending to zero when e — > 0. 

Proof. Thanks to the conservation property for the normalized energy, there exists L{Sq) such that for any 
xq, inf[j„_i,.v„+i] p(x) < L, therefore 

inf u(p)<p)« < C(4). (4.3) 

[xo-R,xo+R\ 

On the other hand, ed^ {eii{p){p)^'^ is uniformly bounded in L\. On the set where p < 1, this is clear. On the 
set where p > 1, by (TC) and Cauchy-Schwarz inequality, 

J <9.,(fi(pKp)^) dx < ^ fi{p)p''-^p,dx 

1 |JeV(p)Wj"^|J„vV'' . C(4,4) 



2,,/„^2|^ 12x1/2 , „ ,1/2 ^, „ ~. (4-4) 



< 
e 



1 The constant also depends on fi, k, but we consider these as fixed. 
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Gathering the two previous inequahties gives iA.l} , from which I I4.2I I follows. □ 

Proposition 4.2 (Higher integrability of the pressure). Fix a compact set K. Under the assumptions recalled at 
the beginning of the present section, for any T > 0, there exists a constant C{K, T, So, Sq) such that 

^ (ppip) + pipK^'ip) + 5k^p))pI) dxdt < C{K, T, So, So). 

Proof. 1. The key commutator identity. Start from the conservation of momentum equation, and integrate 
it over an interval [y, x] to obtain 

I := [pu' + Pip) - ^'{p)u, - (pK^(p)p„ + lipK^'ip) - K^(p))(p.v)')) (y) 

= jA (P") + (P"^ + P(P) - P^(P)".V - ipK'{p)Pxx + 2(Pk"(P) - K'{p)){p,)^)j (X) 

=: II. 

Next, pick a smooth, nonegative, compactly supported function (p equal to 1 on K, multiply the above by 
p{x)cl){x)cf){y), and integrate over y, x e K, and t e [0, T]. We estimate separately the contributions I, II. 

2. Contribution I. Using successively Cauchy-Schwarz inequality, the condition (TC), and Lemma |4A1 it 
can be boimded as follows: 

~ X / P^^^'^^^^'^^^Hp"^ + + P'(P)"? + K'{p)pl + if{p) + pK%p)\p,:^{){y)dxdydt 



<C{T,So)+ r p(x)(|)(, 
Jo 



:^Wy)|<p>"' + P (P)^ ^^\ l\dxdydt<C{K,T,So,So). 



3. Contribution II. This term reads 



II = I ( J\p")) + (P"')(^) + P(PW) - (P'(P)«V) (X) - [pK^(p)p..v - \{pK''{p) - K'{p)){p.f]{x) . 



A3 



Recall that we perform the following manipulation: multiply by p{x)<p{x)(p{y), and integrate over y, x e R, 
and t e [0, T]. Integrating by parts in t, using the equation of conservation of mass, and then integrating by 
parts in x, gives 

^ JJ p{x)(p{x)(p{y)Aidtdxdy 

= Jj p{x)cp{x)cp{y) (pu)dydx - JJ ^(p{x)(p{y)dtp{x) ^{pu) - cp{x)(p{y){pV){x)'^ dtdxdy 



= JJ Pix)(p{x)(p{y) J {pu)dy 

-17/ 



dx 



-d^(p{x)u{x))(p(x)(p{y) 



J ipu) - (piy)(p(x)p{x)^u{x)^'^ dtdxdy 



= 11 Pix)(pix)(piy) ( (pu)dydx + I I I {pu){x)(p'(x)(p{y) I {pu)dtdxdy. 

JJ Jy n Jo JJ Jy 



By Holder's inequality, we find thus 



jj p(x)(Pix)(p(y)Aidydtdx 



<CiK,T,So), 
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while the term gives immediately 

J" J'J'p{x)(p{x)c(){y)A2dydtdx-co J J (ppp{p)dtdx, 

where cg = J (p. Applying successively Cauchy-Schwarz inequality, the condition (TC), and l|4.2|l. one 
obtains 



17/ 



p{x)(p(y)<p{x) A3 dy dtdx 



zio I ^^''^P^^^^'^^'^^'^ ^ J (pp^p)p'^ dtdx 
<C(4) + CC(4),^o) J (ppp{p)dtdx, 
where C is a small constant which will be determined shortly. Finally, a few integrations by parts yield 

J'J'p{x)(p{x)(p{y)A4dxdyds- ^cq J J p{x)pl{x)(^p{x)K'^'{p{x)) + 5K.^(p{x))<p{xfj dx 

+ Co J" J (p'{x)p{x)^K^{x)px{x)dtdx. 

We bound using Cauchy-Schwarz inequality, (TC) and Lemma l4!l1 

r r (p'{x)p{xfK'{p{x))px{x)dtdx < f f \(p'{x)\K%x)p:,{xfdtdx+ f f \cp'{x)\K^p{x))p{x)Utdx 
Jo J Jo J Jo J 

<C(4) + j (l}{x)p'\p{x))p{x)dtdx<C{<^o)- 

Gathering the previous results gives 

X if Udydtdx - '^0 J (pp(p) + ^p(pK"{p) + 5k' (p))plj (x) dtdx 

<CC(4,4)J|" j pp{p)dtdx + C{So,So). 

Taking C to be small enough and combining this inequality with the contributions I and II gives the desired 
result. □ 

4.2 Higher integrability property for the velocity 

General entropy pairs are obtained by integrating the fundamental entropy kernel 



x{p,u;v) = [p^^ - {v - uf\^^ 

with 6 = ^-j- and A = 2(v-\) ^^'^ fundamental entropy flux kernel 

o{p, u;v) = [dv + (1 - 0)u) (p29 -(p- uf 
against an arbitrary function = namely: 



(4.5) 



(4.6) 



if(p,u):= I x{p,u;v)ip{v)dv, q'l'ip,u)-- I a{p,u;v)ip{v) dv. (4.7) 

Jk Jr 



With the choice ip = v\v\, we will simply denote 

Tj-- if, 
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It is convenient in the following to consider that all functions depend either on u, p, or on m, p, depending 
on the context. When derivatives are taken, we adopt the following convention: derivatives in m, u are 
always taken by keeping p constant, while derivatives in p are always taken by keeping u constant. 

The following pointwise bounds hold for and 'q: 



q > p\u\^ + 
\rjm\ < \u\ + p^, 



m < pu^ + p^', 



Ylmp\ < p 



e-i 



\i]mu\ < 1, 



(4.8) 



These bounds appeared first in Lions, Perthame and Tadmor ll49l , and were already used in [47] and next 
in iTlSl . In l24l , these estimates will be checked to hold for a broad class of pressure functions, so that the 
estimate in the following proposition will also be established for general pressure laws. 

Proposition 4.3. [Higher integrability of the velocity] Under the assumptions stated at the beginning of the present 
section, and assuming furthermore that 

ju(p) < p^^^ for p small 
for any T > and every compact set K, there exists a constant C(T, K, (oq, S'q) such that 

J (p\u\^ + p^'^^)dtdx<C{T,K,,ffo,^o)■ 
Proof. The conservation law for the entropy i] reads 

dtrjip, u) + d;q{p, u) = Tj,„{p, u){p'{p)u^)_, + rf,n{p, u) ^pK''(p)p.v.v + ^(pK^'(p) - K^(p))(px) 



HI 



IV 



Multiply this equation by sgn(x - y) and integrate it on (f, x, y) e [0, T] x ]R x K. We examine separately the 
contributions of I, II, III, and IV. 



Contribution I. By I l4.8ll , this term can be bounded by 

I II s^{x - y)\dtdxdy < I I \i]{J,x)\dxdy + I I \j]{Q,x)\dxdy 
I Jo Jr Jk Jk Jr Jk Jr 

< I (pu^ + pA{T,x)dx+ I (pu^ + pA{0,x)dx<2<^o■ 
J^R Jr 

Contribution II. By | |4.8|| , this term contributes 

J I I sgn{x - y) 11 dxdyds = - I I 'q{p{s,y),u{s,y))dyds<- I I (p\u^ + p'^'^A dy ds. 
Jr Jk Jo Jk Jo Jk 

Contribution III. By integration by parts, this term can be written as 

I I I sgn(x-i/)IIIdxdi/ds = - I I I sgn{x - y)(?j„,pp^ + rfnuUx){fi''{p)Ux)dxdyds 
Jo Jk Jk Jq Jr Jk 

+ I I rfm{p,u)p%p)Uxdyds, 
Jo Jk 

which we bound using l l4.8t , and Cauchy-Schwarz inequality: 

J I I sgn(x - y) 111 dx dy ds 
Jr Jk 



(4.9) 



nil 



1112 



1113 



UI4 
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Here, C is a small constant whose value will be fixed in the following. We now bound one by one the terms 
nil to III4. The first one is easy: 



^ J IIIidsdx=^ J ^f{p)\U:,fdsdx<S'o. 
Dealing with the second term is not more difficult, since 

J in2dtdx = j ju'{p)p^^-\A'^dsdx = j ^ ^P^P ''P^ p^i^dsdx < Iq- 
Next, for the third term, we resort to (|4.1| | and the assumption ii{p) < p^l^ for p small: 
^ Uhdtdx = cJ^ J iKi-i^pM^ dtdx 

= (,[ [ lKlu^{p)\ufdtdx + Qf f lKi.i'{p)\ufdtdx 
Jo J\p<l} Jo J\p>l] 

<c r r iK^'{p)p-^'^p^'M^dtdx + [,c(^o,^o) f f iKipy^'M^ dtdx 

Jo J{p<l] Jo J{p>l] 

T 

3, 



< Ce r r p\u\' dtdx + QC{K, T, 4, 4)- 
Jo Jk 

Finally, using once again ||4.1| | gives 



(4.10) 



r f llh dtdx ^ f f 1.1' ip)p^^ dtdx <C{S'o,S'o) f f (p)^ dtdx < C{T,K,S'o,S'o). (4.11) 
Jo J Jo Jk Jo Jk 

Contribution IV. By integrating by parts, this term can be written in the form 

J I I sgn(x - y) IV dsdxdy 
Jk Jk 

^ ~ Jo f f ^^^^ ~ '^^ '^"'""■') (P'^'^(P)P--'-^ + ^(PK^'(P) - K^(p))(p.0^j dsdxdy 

+ Jjjmipis, y), u{s, y)) ^pK^(p)p.v.T + ^(pK^'(p) - K^(p))(px)^) dsdy, 
which we bound using successively (|4.8|l and Cauchy-Schwarz inequality: 

J'J' sgn{x - y)lYdxdyds 

The first four terms on the right-hand side have already been bounded when treating III, and they contribute 
j (iKffipy + fi'{p)ul + lKf(^(p)p''' + f(^(p)p2^-Vx) dsdx < C(T, K, 4, 4) + PM^ ^sdx. 



As for the last two terms, they can be dealt with using (TC): 



I Jo J l^^^^-^ "7^^^^-^-^ do J l^^^^-^ +^^(p.))^Mx 



< ^{<ffo + 4). 



(4.12) 



The desired conclusion follows by gathering the contributions I to IV and taking C small enough. □ 
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4.3 Existence with finite energy solutions to the Euler equations 

Consider solutions of (NSK)^ for which l|1.22|l holds. We gather below all the uniform estimates (uniform 
in e) which have been proved so far (with the subscript e removed for ease of reading): 



1 1 

-pii^ + pe{p) + -K^{p)pl 



^'ip)ul\l, <C(^o,4) 

f^(p)P'(p) 



<C(4,4), 



<C{S'o,S'o), 



(Pxx) + , \Px) 

p p 

y.%p)p 



< C(4,4), 



vip) 



<C(4,^o)o(l), 



||f'(p)<P)1r~ <C(<lb,4)), 

Ha- 

\\pp{p) + p (pK^'ip) + 5k^(p)) (p,)2||^, < C(^o, 5b), 

t,x loc 

\\p\uf + py^% <c(4),4). 

t,x loc 

(4.13) 

We now turn to the proof of Theorem ll.2[ which relies on the energy and higher-order integrability estimates 
stated in | |4.13|| in combination with the compactness framework established in Ii47il for polytropic fluids. 

The entropy-entropy flux pairs {i]^',q'l') were defined in Section l4!2l and we also keep our convention 
on differentiation with respect to u, p, or m, as defined ealier. First of all, we observe that, for any smooth 
and compactly supported function i/', the entropy pair (r]''', q^') satisfies the estimates 



\qHp,u)\<p{p)-', \nl{p,u)\<{p)-', 

riUp,u) <{p)-', nip{p,u) <p'-\p)-', 



(4.14) 



as is easily checked from the expressions ( |4.5l l-l [4l6l l of Xr o^- Iri order to establish the convergence property 
stated in Theorem 1 1.21 it now suffices to establish the following result. 

Proposition 4.4 (Compactness of the entropy dissipation). Assume that the tame condition (TC) and the 
bounds i4.13[ hold. Then, for all compactly supported test-functions \p, the sequence dtri^{p'^, if) + d^q^ip^, u^) is 
compact in H~^J^+ x ]R). 

Proof. 1. Conservation law for f]'''. It reads 

diTf'ip, U) + d:,q'i'{p, U) = rji (^{p'{p)u^)^ + (pK^(p)px.r + ^ipK^' (p) - K'{p)){p,y 

We fix from now on a compact set K, a time T > 0, and we wiU show the compactness of dtif'{p^,u^) + 
d^q'i'ip', u') in H-i([0, T] X K). 



2. Compactness of the viscous term in W^^^''' for q < 2 We deal first with the viscous term r]1j,{p'^{p)Ux)x, 
which can be written 

rit{p^{p)u,)^ = (r]ip'{p)U:, - r]ippxp''{p)ux - rituUxP%p)u^ ■ 
I II in 

Using successively the estimates (|4.14|l and the bound on ^-^^ £ i^J",./ we obtain 

l^'{p)Uxr]t ^ P'ip) l".vl $ ^|p'ip)\ux\ Vf'^(p) = 0(1) ^ffFip)\uA 

where o(l) is understood in the L^^. topology and e ^ 0; namely, p^{p) = 0(1) follows from p^{p) < 
inf(e/j(p), (p)~^). Therefore, we find 

1^0 in H-H[0, T] X K)) as n ^ +co. (4.15) 



J', 
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(4.16) 
(4.17) 



Next, in view of the estimates 14. 14l and by Cauchy-Schwarz inequality, we have 

|ii| < p"-V'(p)|w..llp.vl < ^'{p)\uA^ + p^'-^ii'{p)\pf < ^Hp)\uA^ + ^^-^^^^pl 

This implies a uniform bound for II in the norm: 

l|II|lLi([0,r]xK) ^ 1- 

Finally, by Lemma 14.141 we have |III| < iu'^{p)ul, thus 

I|III|Ili([o,t]xk) ^ 1- 

Gathering (|4.15l l, (|4.16l l, and (|4.17) l gives that rit,{n(p)ux)x is compact in W"^'"? for all q < 2. 

3. Compactness of the capillary term in W~^'i for cj <2 We now deal with the capillary term 7]fj{pK'^(p)pxx + 
|(pK'^'(p) - K^{p)){px)^)x which can also be written as 

vfn{pK^ip)pxx + ^ipK^'ip) - K^p)){pxf)x = {ntpK%p)pxx ^ ' pK^p)pxxr]fnpPx - pK^ {p)pxxT]fnuUx 

i 11 III 

+ I {ntipK^'ip) - K^{p)){pxf - UpK^'ip) - K'{p)){pxfrit,pPx - UpK^'ip) - K^{p)){pxfl]tuUx 



IV v 
Using the estimates l|4.14|l and the condition (TC) gives 



VI 



riipK''{p)pxx < pK'pxx < yU^Hp)-)J—^pxx = o(l). 



Similarly, the estimates (|4.141 , the condition (TC), and Cauchy-Schwarz inequality yield 

|II| = \pK%p)pxxntppx\ < pK'{p)\pxx\p''-'\px\ < ^{Pxxf + ^(Pxf- 

Next, the estimates (|4.14l l, the condition (TC), and Cauchy-Schwarz inequality give us 



mil = 



pK^{p)pxxr]tmUx < \pK'{p)pxxUx 



{pxxf + P-'iUxf, 



To deal with IV, we use (|4.14|l and the condition (TC) to obtain: 



r]l{pK''{p)-K'{p)){pxf <K'{p){pxf< ^[i^)J^{pxf = o{l) 



It^iPxf. 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



Once again, (|4.14b , the condition (TC), and Cauchy-Schwarz inequality give 
|V| = 

Finally, the same arguments give 



(pK^'(p) - K^p)){pxrr^lppx < iK^p'^-'ipxYpxl < -^{Pxf + ^{Pxf, 



P' 



|VI| 



ipx^'ip) - K'{p)){pxfntuUx\ < \k'pIux\ < lf{Uxf + "^{Pxf. 



(4.22) 



(4.23) 



The estimates l|4.18|l and l|4.21|l imply, together with the estimates listed in l|4.13|l , 

I + IV^O rnH-i(Kx [0,T]))asM -> -l-oo, 
whereas the estimates l|4.19ll , l|4.20t , l|4.22t , and l l4.23ll entail, together with the estimates listed in l|4.131 , 

||II + III + V + VI||ti(Kx[0,T]) 5 1- 

Combining the last two statements implies that the capillary term is compact in W"^''' for q <2. 

4. Interpolation argument and conclusion. In view of l|4.14|l , the functions |r]'''| and \q^\ are bounded by 
p. Since p is boimded in U^^, it follows that dtif^' + dxq^ is bounded in W"^'' "^^, with of course y + 1 > 2. 
Interpolating this property with the compactness in W"^-' for aWq <2 gives the desired result. □ 
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4.4 Completion of the proof of Theorem 11.21 



Following 113, we can associate to the sequence (p"^, m*^) a Young measure r' : ]R+ x ]R ^ Prob(]R+ x M), 
satisfying the weak convergence property l ll.31ll for all continuous functions / = f{p, u) satisfying a certain 
growth condition. More precisely, the higher-integrability property of the density allows us to check that 
jl.Sll holds for all |/| < /o(l + p''"^^) with limp^+cw /o(p) = 0. This is sufficient to imply that the entropy 
rj'l'ip^, u^) and the entropy flux q'^'ip^, u^) converge to {v, if') and {v, q^') , respectively, provided the function 
ip is compactly supported. This convergence property suffices to state Proposition l4.4l above. 

However, in order to recover the conservation laws in the Euler system and the global energy inequality, 
the additional integrability property for the velocity is required, which allows us now to use sub-cubic 
functions ip for dealing with the entropies, and sub-quadratic functions ip in the entropy flux. Observe that 
the local energy identity does not make sense at the level of regularity and integrability imder consideration 
in the present paper. 

Finally, by applying the reduction lemma established in [47) for Young measures satisfying Tartar's 
equation and associated with polytropic fluids, we conclude that v is a Dirac mass or else is supported on 
the vacuum line. This completes the proof of Theorem ll.2l We refer to | [47| as well as |25] for further details 
and generalizations, including a framework covering real fluids. 
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